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ABSTRACT 


This study undertakes to improve the theory and numerical analysis of the ancient 
astronomical observations (-1374 to 1715), and to combine these with the modem 
data in a simultaneous solution for: J’it: the tidal acceleration of the lunar longi- 
tude; (I)/a) the observed apparent acceleration of the earth’s rotation; the 

true nontidal geophysical part of tli/w; and 6/G the rate of change in the Gravi"’ 
tational Constant. 

Error in the lunar node is shown to have a large systematic effect on ftom 
eclipses, and corrections of +4.39 "cy”^ to and -0.10 "cy“^ to Q are added to 

the standard epheme rides based on the results of other studies. The complete 
formulation for AT = A + B’T + C*T^ is essential to finding d/o), and in testing 
the constancy of tli/o) and during the historical period. An exact linear 
inequalities Statistical processing of the large solar eclipses is substituted 
for traditional least squares in order to preserve observational sensitivity in 
the separation of and A strictly observational separation of and 

G/G is made by including ancient data on Universal Time (UT) and very modern data 
on Atomic Time (AT) in the solution. 

The ancient eclipses alone provide « -30.0 ± 3.0 "cy“^ and the corresponding 
d)/aj = -24.3 ± 2.0 X 10"^^ yr*”^ , and AT = 20 + 114-T + 38. 30*T^ (T in cy from 1900) . 
This (L/oj agrees very well with dj/w from timed solar equinoxes alone (independent 
of the moon): oj/lo - -23.8 ± 2.3 X l0~^-yr“V and gives an independent verification 
that the from eclipses is not significantly biased. The observational bounds 
imposed on AT as a function of time prove that the accelerations have remained 
constant during the historical period within uncertainties. These three results 
from the ancient observations are then combined with: (i) = -26 ± 2 "cy”^ f rom 

transits of Mercury (1677 - 1975) and (li) = -36 ± 5 ”cy“^ from lunar occulta- 
tions on Atomic Time , provided by other author's studies. The simultaneous solution 
from all data provides ~ -27,2 ±1.7 "cy~^, di/O) = -22.6 ± 1.1 X 10“^^yr"^» 
AT 5= 20 + 114*T + 35.55*T^; where, on the so-called primitive cosmologies 0?.^* 
Hoyle-Narlikar) , G/G = -2.3 ± 1.5 X 10“^^ yr“^ with between +3.1 and +5.9 

± ^^4 X 10“^Vyr~^ depending upon the assumed earth expansion f-actor as a function 
of S/G; or on the Dirac cosmology, &/G = -5.1 ± 3.0 X 10“ ^^yr-^ and = +3.1 

± 4.4 X 10”^^ yr'“^ ; or constraining G/G = 0, - +9.2 ± 2.5 X I0“^^yr“^. The 

Brans-Dicke cosmology is consistent with a "primitive" G/G as found above (near 
-2. 3 X 10'^^ yr’*^). In the remaining cosmologies, a Hubble Constant of 55 ± 7 
km/s/Mpc corresponds to G/G = -5.6 ± 0. 7 X 10“^^ yr” ^ The Dirac cosmology agrees 
well with this result from extragalactic data, the others needing more (h/a>pj^X 
reveal G/G than is available. VanFlandem finds G/G from data sets (i) and (ii) 
above: G/G - -5 , 8 ± 3. 0 X 10"^^ yr“^. Solving for G/G on any cosmology above gives 
^/^NT consistent x^yith zero. Seve^xal possible (but unproven) geophysical theories 
have been advanced to explain a in the range found above, including post 

glacial uplift, geomagnetic torques, and phase- change core. Assuming uj/a>j>{gX = 0 
with the timed ancient equinoxes' (b/w above, independent of (ii) and the moon, and 
data set (i) , provides an estimate of G/G ^ -6.9 ± 3.0 X 10“ yr“^ independent of 
VanFlandem' s Atomic Time data and consistent with the other determinations. It 
appears that either we really have a cosmological G/G consistent with the Hubb le 
Constant, or we have a significant (b/ajj-j^-j. The former is consistent with all 
observations and the relationships between 6/G and the Hubble Constant contained 
in several modem cosmologies; the latter can be forced only if the data (ii) is 
wrong we abandon the Hubble Constant versus G/0 relationship common to modern 
cosmologies . Final with Hubble rate on Dirac is: +2.5 ± 2. 5 X 10“^^ yr“^ . 

This provides three independent:^ determinations of a rate of change of G consia- 
ten t xs^i th the Hubb le Cons taitt and a near-ze ro non t idal rot ation al acce Te vat ion 
of the earth, The tidal accelerations are shoxm to have remained constant 

during the his torical period within uncertain ties . Ancient an d modem solar system 
data, and extragalactic observations provide a completely consistenv: astronomical 

and cosmological scheme. 
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Figure 1 


This paper is an updated summary of the prin- 
cipal i-esults of Muller (1975)^ The study emphasizes 
the theory and numerical analysis of the ancient 
astronomical observations, taking the historical and 
linguistic background from other authors. The 
strongest data are the large solar eclipses. The 
lunar parallax (distance) was determined from such 
an observation by Hipparchus (oa. 135BC) . Sir Edmund 
Halley (1695) was the first modern analyst to use 
ancient eclipses, and he discovered the so-called 
secular acceleration of the moon’s longitude. His 
discovery was confirmed and refined by Dunthome 
(1749), Baily (1811), Airy (1857) and Newcomb (1875) 
among others before the present century. 

Euler (1770) won the prize of the French 
Academy for a theorGitioal explanation of the secular 
acceleration with a proof that it could not come 
from any defect in lunar theory. Lagrange (1774) 
proved that the figure of the earth, and hence the 
planetary terms, could not be responsible. Laplace 
(1786) found a large theoretical contribution from 
the changing eccentricity of the earth’s orbit, 
proving that Euler had erred, but this accounted for 
only part of the discrepancy. The modern lunar 
theory was created by Hansen (1857) and Brown (1919). 

All this time it was assumed that the secular 
acceleration of the moon was veal as opposed to 
appavent, Spencer Jones (1939) proved that part of 
the acceleration Was apparent only, arising from the 
slowing of the earth’s rotation. The lunar longi- 
tude, and earth’s rotation, were both accelerated by 
the effects of tidal action as suggested by Immanuel 
Kant (1754) nearly two centuries earlier! Kant made 
the astounding argument that the earth’s rotation 
rate could not be constant because the tidal drag 
must inevitably slow it down! Kelvin (1897) notes 
that nobody anticipated Kant in this fundamental 
insight, and Hastie (1900) provides a translation of 
Kant's treatise with a valuable commentary- 
Amazingly, Kant later (in his Phyoical Geography) 
withdrew the argument, and suggested that this 
effect might be offset by the accretion of the 
earth’s core! Urey (1952) considers a similar 
hypothesis , and it is still being debated by geophys- 
icists , along with other possible sources of change 
in the earth’s moment of inertia (as below). 

The so-called Spencer Jones anomaly emerged as 
modern studies yielded values of the lunar tidal 
acceleration from ancient observations (top of 
Figure 1) and Wdem observations (second block of 
Figure 1), which differed substantially. Newton 
(1970) argued from his data analysis that the lunar 
tidal acceleration had changed during historical 
time, but this view conflicts with theory as noted 
by Hunk & MacDonald (I960) and with the data analysis 
in this paper. Both the lunar longitude and earth’s 
rotation are accelerated by the tidal couple, and 
both can be independently determined from the 
observations in the manner described below. Since 
angular momentum is conserved, we can check the 


Previous Determinations of Lunar AcGeleration 


Author 

Epoch 

h 


an 

Fotheringham (1920) 

-200 

-30. 

.8*^ 


De Sitter (1927) 


-37. 

7 

4. 

3 

Newton (1970) All Ancient Data 

-200 

-^1, 

6 

4. 

3 

Newton (1970) All Medieval Data 

1000 

-42. 

3 

6. 

1 

Stephenson (1972) 

-300 

-34, 

2 

1. 

9 

Muller & Stephenson (1975) ILE node 

-400 

-37. 

5 

5. 

0 

Muller (1975) Ancient eclipses, ILE 

-470 

-34. 

5 

3. 

0 

This Paper: Eclipses, corrected node 

-470 

o 

cn 

1 

0 

3. 


Spencer Jones (1939) 

1800 

-22. 

4 

1. 

1 

Clemence (1948) 


-17. 

9 

4. 

3 

VanFlandern (1976) 

1965 

-35, 

0 

5. 

0§ 

Oesterwinter & Cohen (1972) 

1940 

-38. 

0 

8. 


Newton (1968) Satellite Data 

1965 

-20. 

1 

2. 

6 

Morrison & Ward (1975) 

1830 

-26, 

_0_ 

h 


Lambeck (1975) Prom Tidal Theory 


-35. 

.0 

4. 

0 

"cy-2 

Newton (1970) finds the corrected value -34.0 
^ On Atomic Time and will include effects of G/G 
VanFlandern (private communication) prefers range • 
to -38” cy"^ ; some of the data are on Atomic Time. 

-18 


observational determinations for consistency. 

Results of previous studies, Newton (1970), 

indicated a very large residual nontidal acceler- 
ation in the earth’s rotation d)/u)jq^ - 23 X 10" ^ Vr“ - 
This must arise from changes in the moment of iner- 
tia or interchange of angular momentum inside the 
earth, or a rate of change in the Gravitational 
Constant as first suggested by Dicke (1957) i The 
present research was primarily motivated by the 
existence of these two long-standing unexplained 
anomalies . 

Theoretical Foundation 

Ne are engaged here in a comparison of lunar 
(and Solar) orbital theory with observations of the 
sun and moon. The reference lunar ephemeris is the 
Improved Lunar Ephemeris (1954), often denoted "ILE”. 
The theory of the sun is that of Newcomb (1895) . 
Corrections derived in this study are to be added to 
the mean elements , nominal values of which may also 
be found in the Explanatory Supplement to the 
Astronomical Ephemeris (1961) . The theories are 
constructed in such a way that if you desire a 
position, you begin with the mean values, and then 


add the theoretical periodic perturbations. The 
expression for the mean lunar longitude is: 

Xm = Ao + Bq-T + (Co + Ct-)-T'^ + D^-tS (1) 

where T is expressed in Julian centuries from 1900, 

0 January, 12 hours ET (Ephemeris Time)* Eor the 
lunar node: 

n = Eq + Fq-T + +Ht*T^ (2) 

Eor the mean solar longitude Xgi 

= lo + + Kt-TV+ Lt.T^ (3) 

The letters A through L represent the ephem- 
erls coefficients for the mean positions at epoch 
and the subscripts *'o" and ‘*t” denote observation- 
ally and theoretically determined values respect- 
ively. It is clear that A, B, E, I and J must be 
determined by observation since these are the 
positions and rates of the parameters at the epoch 
(which is quite arbitrary). Tlie laws of Newton 
relate B and J to quantities such as masses of the 
primary bodies and the distances of the orbiting 
bodies, but these must be observed. If there were 
no tides acting on the earth-moon system, then the 
coefficient C would be strictly determined theoreti- 
cally from Newton’s theory. The tidal part must be 
observed as an acceleration in lunar longitude, Gq, 
and adds to the theoretically computed part 
That determination is central to this analysis and 
many before it (figure 1). The earth’s angular 
momentum about the sun is so large that tidal 
effects on coefficient K are negligible. 

Since "n'Ms the standard symbol for the lunar 
(or planetary) rate in longitude, n Is the acceler- 
ation which we must observe, Subscripted ”t" for the 
true tidal part, and "a" if observed on Atomic Time 
instead of Ephemeris Time, where: 

^o ^ 

As first noted in Di eke (1957) , these two time 
scales will differ if the Gravitational Constant or 
an e q ui valen t phy s i c al p ar ame te r i s ch an g in g , r ath er 
than remaining constant as assumed in Newtonian 
mechanics. This is considered further below. 

Problem of the Lunar Node 

In principle, (2) should be theoretically 
calculable from Newtonian mechanics , given the 
masses, orbits and figures of the earthy moon, Sun 
and major planets. As noted in Eckert (1965) and 
elsewhere, this has not been a practical undertaking. 
An observational determination due to Brown (1914) 
and reconsidered by Spencer Jones (1932) was used 
in the ILE (1954). Muller, Newhall, Van Flan dern & 
Williams (1976) have very recently run a computer 
processed numerical integration of the solar system 
Including the moon back more than 250 years. 


Although initial orbital conditions and physical 
constants such as masses and distances had to be 
observationally determined in other research, 
primarily lunar, planetary and spacecraft laser and 
radio ranging, it is only necessary to have a 
reasonably converged integration to determine the 
lunar nodal rate. The result is given here as the 
correction to be added to the nominal ILE (1954) 
coefficient F^: 

f2 = + 4.39 + 0.15 "cy"^ (5) 

The formal uncertainty is provisional and is still 
being studied as work on the referenced paper is 
completed. This result agrees very well with the 
observational value of Martin & VanFlandem (1970) 
of ft-+4.31±0.36 ”cy~V, which VanFlandern (private 
communication) had worried might be adversely 
affected by systematic errors in the observations; 
and with Muller (1975) from a combination of ancient 
and modem observations similar to that undertaken 
below, of = +7. 1 ± 3. 4 “cy”^ (assuming &/G-0), and 
+ 8.2± 3,6 ’’cy*”^ (solving for G/G). The correction 
in (5) will therefore be made to the ephemeris (or 
results) used in the reduction of the ancient 
eclipses . 

Lunar Node as a Systematic Error Source 

Newton (1970), page 287, dismissed node as a 
])oss±hle systematic error in the analysis of ancient 
observations on the grounds that the small effects 
(5) average zero. His argument is sound for the 
timed observations, but Muller (1975) showed that it 
is incorrect for un timed data such as large solar 
eclipv'es (the dominant data source). Cowell (1905) 
had suggested that node errors would not be separ- 
able from the lunar longitude acceleration n in 

Figure 2 

Equivalence of Node & Longitude Change at Conjunction 
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ancient observations, but his point was overlooked 
in the studies, figure 1, before Muller (1975). 

If the lunar orbit was exactly in the plane of 
the ecliptic, the eclipse observations would-be 
unable to separate y\ and J)/w (the earth’s rotational 
acceleration). This can be seen in Figure 2 by 
noting: [i] u)/w cannot be distinguished from an 
acceleration of the sun’s longitude; [ii] neglecting 
lunar parallax change and other very small second- 
order effects arising from changing the sun’s 
longitude, the acceleration of the moon cannot be 
distinguished from that of the sun if their orbits 
are in the same plane* The lunar inclination I of 
about 5® provides latitude variation A^) at conjunc- 
tion which maps to the earth’s surface in geometric 
eclipse limits (for example). The altered geometry 
provides the small, but essential, separation of A 
and d)/U). 

For a nominal conjunction at the node Cq in 
figure 2, a change in A which would move the moon to 
the left of node an amount AX, as shown, gives rise 
to a new conjunction to the right of node at due 
to the sun’s mean motion ns* This happens because 
the moon will reach the node sooner, and the sun 
will not yet have reached it. The sun is at Sfi and 
the moon is in conjunction at and : ‘ 


became available only a few days before this 
writing, it will be simplest to use the calcula- 
tions of Muller (1975) and apply a correction to h 
based on The change in apparent resulting 
from a change in epheraeris nodal acceleration AQ 
follows from (8): 

AA ^ -12.15AS (9) 

For an equal Af2 from rate and acceleration terras we 
have Combining this with (9) yields: 

M = {2(12.15)/T}^! = 24.3«/T (10) 

This gives the effect on solved from eclipses at 
mean epoch T in the presence of an error in nodal 
rate ft, The ancient eclipses in Muller (1975) which 
participate in the mean acceleration solutions range 
from the year 120 back to -1374, figure 14 below, 
and have a weighted mean epoch near the year -470. 
This epoch can be used with sufficint accuracy in 
the Correction to n arising from ft, and substituting 
T=23.7 into (10) we have the change in apparent An 
resulting from a change in ephemeris nodal rate Aft: 

An = l,02cy"*AQ (11) 


A4) = ± AX ( 


nj^CosI 


-} Tan I 


( 6 ) 


A direct motion of the node itself gives rise to a 
conjunction at Cq with the sun at sft and the moon at 
Mft, and: 


A<!) == ± AQ Tan I 


(7) 


where the (-) sign applies at the ascending node and 
the (+) sign at the descending node. The key point 
is that in untimed observations we can ignore the 
time difference implied by the sun’s position at 

or and changes arising from the incTination of 
the earth's equator, and moving the node and longi- 
tude aZwaijB gives rise to the same geometric effect 
A^ at each node. Furthermore, because of the sun’s 
relatively small mean motion ng compared with the 
moon’s a change in the node is about 12 times as 

effective in changing A4> as is a change in longitude 
(Aa-AX): 

Acj) (from AJ2) *= 12.15 A^) (from AX) (8) 

There are other second-order variations Which affect 
the projection of the A(j) at conjunction to the 
earth’s surface in eclipses, but they are about a 
factor of five smaller than the basic geometric 
term in (8). A numerical test of this analysis can 
be easily made by inspecting the ratios of the 
partials with respect to A. and ft for each observa^ 
tion in Figure 9 below. 

Muller (1975) caTculated eclipses assuming 
ft - 0, but anticipated that there would probably need 
to be a correction when a reliable f2 became avail- 
able. Since the apparently definitive result (5) . 


or for the correction implied by (5) : 


An ■t 4.5"cy“^ (12) 

which can be added to the uncorrected h found on the 
standard ephemeris as necessary below. 

Muller (1975) undertook a reexamination of the 
basis for ft in the analytical ephemeris. The 
coefficient in (2) x^as computed by Brovm (1919) . 
Deprit & Henrard (1975) provide the partials 
relating various lunar parameters , which VanFlandem 
(private communication) has verified from the 
classical analysis of Delaunay. If the planetary or 
lunar orbits are (for example) as 

revealed by the apparent longitude accelerations, 
there is an effect on the lunar nodal rate: 

9ft 9ft . 

— =+.00375; — = -.1038+ (13) 

3ns 

wh i ch yie Ids the co r re c t i on t o n o d a 1 ac ce le ration : 

jj = ^ + ^ = .00375Aa - .lOSSAg (14) 

I f th e orb i t s are exp an ding , then th e n od a Irate 
must change accordingly , since it arises from gravi- 
tational perturbations. This means that Brown (1919) 
has neglected the effect of the tidally induced 
expansion of the lunar orbit on the node. For an 
near -30" cy”^ , equations (14) and (9) yield an 
additive correction of about +1.5 "cy”^ to the 
apparent n^. on the standard ephemeris. Tliis 
correction to ft (and hence nt) was made in Muller 
(1975) and is therefore already re fleet e d in the 
Misprinted +.1559 in some Deprit & Henrard preprints. 


tabular computations of the eclipses below. When 
the final value of is adopted in the lunar 
ephemeris , the equivalent change in required by 
(14) must be entered in the coefficient G(. of nodal 
acceleration (2). 

The question of a rate of change of the 
Gravitational Constant is subtle when it comes to 
effects on the lunar node. The partials in (14) 
from sun and moon are such that they nearly cancel 
when both orbits expand proportionally, and nodal 
effects are rather small. This question will be 
considered below, with suitable coefficients placed 
in the equations of condition. 

The Earth *s Rotational Acceleration 

Since the earth rotates on its axis faster than 
the moon orbits about the earth, the nearer tidal 
bulge on the earth is dragged ahead of the moon. 

The torque on the moon arises from the net differ- 
ence between the gravitational attractions exerted 
by the moon on the two tidal bulges. Figure 3. 

Figure 3 

Earth ~ Moon Tidal Couple 



The nearer bulge is larger* and closer, and the moon 
is accelerated (gains energy), and moves farther 
from the earth, slows down, and exhibits a n&gdtiVe 
apparent acceleration in longitude. The earth’s 
rotation is also negatively accelerated. 

An exceptionally clear analysis of the differ- 
ent time scales used in Astronomy is to be found in 
Mulholland (1972). Ephemeris Time (ET) is the 
independent ‘argument of the ephemeris required for 
this analysis - equations (1) through (3). Universal 
Time (UT) is the measure of the rotational position 
of the earth with respect to the sun defined by 
equation (3). The earth’s rotational behavior is 
measured by examining the function AT defined by: 

AT ^ BT - UT (15) 

as a function of time (e. 0. from -1373 co 4''1975) , 
where ET and UT are the readings on the Ephemeris 
and Universal Time "r.locks" at Various instantaneous 
epochs of observation, 1/ the ephemeris of tba sun 
(3) is correct, then the observed values of AT in 


the presence of an acceleration in the earth’s 
rotation (b/to will satisfy: 

AT = C*T^ where d)/ui * kC (16) 

If C and cj/o3 are in the same units, 2, It has 
been customary to express n in arcseconds per 
century squared, *'cy“^, w/w and G/G in fractional 
parts per year, 10"^^ yr’’^ For simplicity in this 
paper, the expressions for units conversions will be 
given without the awkward units, for example, the 
”k" of equation (16): 

u}/ai=-.635 C (17) 

Previous analyses of the mtoient observations 
assumed that (3) was sufficiently correct, and used 
(16) to define and find oj/o) from AT values at many 
epochs of observation throughout history; see for 
example Newton (1970) page 6, Fotheringham (1920), 
Curott (1966), Dicke (1966) and (1969), and 
Stephenson (1972). Studies of modem observations, 
including Spencer Jones (1939), Clemence (1948), 
Brouwer (1952) and Van der Waerden (1961) all made 
use of the expression defining AT: 

AT = A + B-T + C'T^ (18) 

The A and B coefficients in (18) are equivalent to 
the Iq and Jq coefficients in (3)* and A&B will be 
zero if and only if 1^ & Jq are correct. Muller 
(1975) citing analysis by Martin (1972) shows that 
the data available to Newcomb (1895) were insuffi- 
cient to adequately determine Iq and Jq in (3), and 
the substantial value of coefficient B found here, 
equation (62) below, and in Brouwer (1952) for 
example, is real and essential to finding an 
accurate value of d)/(D from the ancient observations. 
One could correct (3) , but it is customary to leave 
the solar ephemeris fixed, and correct, instead, the 
form (18), and that is the approach used here. 

Muller (1975) and Muller & Stephenson (1975) show 
that the a dop t i on o f the simp li f yin g ass ump t ion ( 16 ) 
by Newton (1970) and (1972) gave rise to his 
apparent observational conclusion that the accelera- 
tions had varied substantially over historical time, 
contrary to the findings here and in the other 
references. This modeling simplification also 
tended to substantially overestimate the true, 
nontidal portion of uj/oj, in Newton (1970) , 

and this parameter will be considered in the next 
section. 

Earth - Moon Ti dal Interaction 

Gonservation of angular momentum in the earth- 
moon system implies that the tidally induced 
acceie rations of the earth’s rotation uj/a)^ and the 
moon’s longitude are related: 

0.935^t (lunar tide alone) (19) 


- 4 - 
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The solar tide must be added to the lunar contri- 
bution. Jeffreys (1962) gives equations for 
computing the ratio between the lunar and solar 
tide, and Newton (1968) evaluated them with the 
results: 5:1 for linear tides, 3.8:1 for quadratic 
(shallow sea) tides. These analyses used approxi- 
mations to the tidal model, and Lambeck (1975) is 
the most recent work and carefv^lly distinguishes the 
contributions of sun and moon, i^aula (private 
communication) notes that Lambeck* a analysis implies: 

1.224^^- (lunar + solar tides) (20) 

Transfer of angular momentum from the solar tide to 
the earth Vs orbit is negligible. There is a signif- 
icant solar atmospheric tide to be accounted for, 
which Munk & MacDonald (1960) estimated at +1. 6 X 
10~^^yr"^. Newton (1970) used +2.7, but Lambeck 
(1975) finds from Kertz et.al, (1970): 

6o)/ti%= + 1.1 X 10”^^yr“^ (21) 

Lambeck notes the error in Newton, and adoption of 
(21) is undoubtedly the best available result. To 
compute the appavent non tidal acceleration from the 
observed total acceleration: 

"" ^ 1 ^ - 1.224j^ - 1.1 (22) 

where the units are those adopted above. 

Observational Determination of the Accelerations 

To extract the several highly desirable astron- 
omical and geophysical parameters from the observa- 
tions, it is necessary to independently determine 
{i.e* separate) the parameters (h/to and n. As Newton 
(1970), page 3, and others have noted, the determina- 
tion of one parameter given the other is done to a 
high precision, whereas the solution as independent 
parameters is difficult and the sensitivity is not 
all that one might desire (even two parameter 
systems can be highly correlated! ) . The solutions 
given in this paper, Muller (1975) , Muller & 
Stephenson (1975), and Newton (1970) page 2*?2 all 
agree in the parameter AT (at a n) within 

about 200 seconds of time (see figure 13 below) , 
lor oj/o) as an independent parameter, the -typical 
uncertainty for an ancient epoch (near 0) will be 
at least an order of magnitude greater in absolute 
At. The high degree of correlation between n and 
(b/m can be seen in figure 2. Consider the effect of 
a change in n of 5”cy”^ at the epoch -100 (T = -20 cy). 
This changes the computed lunar longitude AX at the 
seme ET by ab ou t 10 0 0 ■ ’ . As we s aw ab o ve , thi s can 
be. cancelled by an appropriate change in ET-UT, 
leaving only the small change in latitude at 
conjunction A(j) - 8" from (6). In a solar eclipse 
path on earth, this 8" is equivalent to about 15km 
in latitude change, which for an inclination of 30° 
between the eclipse path limit... and a line of equal 
latitude gives about 30km of longitude change with 


which to measure the Independent AT change between 
the parameters. Total eclipse paths are typically 
100 to 200km wide, so it is apparent that a few 
sound observations can probably resolve to 5 ’’cy""^ 
or better. 

In timed lunar observations, the situation is 
much worse. The earth will rotate the above 30km in 
about 1.5 minutes, while the moon will cover the 
roughly 80'* of residual longitude change at conjunc- 
tion (Cp to M^ in figure 2) in about 2.5 minutes. 

The net effect of these timing variations in a given 
timed observation will be some linear combination of 
these numbers with coefficients of less than unity 
applied to each, depending upon the circumstances of 
the observation, or something in the range of 1 to 3 
minutes typically. I know of only one observation 
from the ancient world which even approaches this 
level of accuracy: the Babylonian timing of the 
beginning of a partial eclipse of the sun on 26 
September -321 as reported by Fotheringham (1935). 

It turns out to have a potential observational 
accuracy of about 5 minutes of time, and a sensi- 
tivity of about 1 minute of time per 4 "cy“^ in ^ 
(3AT/9n in figure 9), It can therefore bound h 
within roughly ±20"cy“^, and it is by far the most 
precise timed observation we have from that era! 
Muller (1975) goes on to justify the conclusion that 
ancient timed lunar observations and untimed lunar 
eclipse data have insufficient sensitivity to have 
any material effect on the parameter solutions. 

These data are therefore not considered further. 

The timed soZar observations (equinoxes) are 
quite another matter since these data are indepen- 
dent of the moon and measure AT directly. They 
therefore separate the parameters with the full 
accuracy of the observations themselves , and a 
precision of an hour can be very useful! These 
data are considered below, in Muller (1975) and 
Newton (1970) . 

The mean linear relationship between the 
parameters can be computed from the mean motions of 
the moon and earthVs rotation, and is a useful 
formula to note. For observations of the moon with 
respect to the sun; 

= +.622; m/u - .622 A = L (23) 

where L is the observed (solved^- for) constant term 
in the linear relationship between the parameters.. 
Individual observations will have partial deriva- 
tives which can differ from .622 by up to about 1052, 
the limit being imposed by the factor Tan I in (6). 
This is another way of viewing the separation of the 
parameters. 

S te ph ens on (19 72 ) , and Mull e r & S t ephe n s on (1975) 
argued convincingly that many partial eclipse data 
were subject to very severe systematic bias due to 
population distributions. To be safe from this 
avoidable bias one must use total eclipses only. 
Muller (1975) also showed that the partial eclipse 
dat a had ins i gn i f i cant sens it ivi ty in sep arat ing the 
parameters, and could safely be deleted. 
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Minimum Deletion Linear Inequalities Filter 


M import an c finding of Muller (1975) was that 
least squares is not the optimum statistical filter 
in the reduction of total solar eclipse data. Every 
study of ancient eclipses, except Fotheringham 
(1920), converted each eclipse observation into a 
weighted equation of condition suitable for least 
squares processing. This Will work reasonably well 
if the desired (or achievable) final solution 
uncertainties are large compared with typical path 
widths. Stephenson (1972) who relied heavily on 
the most ancient data in his mean solutions, and 
Curott (1966) who sought limited accuracy, lost 
comparatively little by this statistical simplifi- 
cation. There are two imneecasapy losses in 
sensitivity which result from using least squares 
equations of condition with total solar eclipses. 

[1] The hard limits (see below) of the path edge are 
approximated by the Gaussian distribution Curve and 
lose sharpness. [2] The estimator tends to force 
solutions towards the center of paths, when stations 
of observation within the path are equally likely. 
Muller (1975) found Newton’s (1970) solutions to be 
capable of considerable variation: ”lt has been 
possible with Newton’s data to obtain reasonably 
stable solutions for n as an independent parameter 
in the range -29 to -46 "cy“^.’' Newton (1970) and 
(1972) largely dissipated his accuracy through this 
approach because he placed heavy demands on the 
solutions including higher precision and an attempt 
to see if the accelerations varied (necessitating ¥ 
solutions at several epochs) . Muller (1975) and 
this paper use the path limits (not just widths) 
as bounds on the possible variation of the param- 
eters. Such a situation is unusual in observational 
science, but not unknown, and the exact mathematical 
approach to the statistical reduction of such data 
is to use systems of linear inequalities. 

Figure 4 

Observational Equations of Condition 


N 



Figure 4 illustrates the parameterization of a 
total eclipse observation. The eclipse path is 
calculated for the nominal parameter values. 

Changing (h/a) moves the computed path directly east- 
west, whereas changing A moves the path in a 
direction which is generally not east-west, differ- 
ing by the A(|) at conjunction noted in figure 2. If 
an observation places station “C" within the path of 
a knom eclipse, both upper and lower limits are 
placed on the maximum range of the two parameters, 

Tf N and S are the values for the clvange in the 
parameter which brings the observation station 
onto the north and south limits respectively, 
and for 6n, then figure 4 gives the inequality 
equations of condition on the parameters. There 
are corresponding equations for an observation 
requiring the station to be outside the path (strict 
denial of totality). Each parameter is bounded by 
two equations of condition, which yield a pair of 
linear hOimdarii^c corresponding to each equation in 
the solution space for the two parameters, which in 
principle can be plotted as in Figure 5. 

Each equation of condition implies that values 
of the parameters on one side of the line (cross- 
hatched) are ruled out. The logical intersection of 
these constraint equations is defined, as the 
solution space, within which the true solution must 
lie (assuming correct observations). Figure 5 is 
similar to the presentation of Fotheringham (1920) , 
although he did not use an orthogonal pair of 
parameters. 

A given set of linear inequalities will either 
provide a solution space with finite area, or the 
set of equations will be inconsistent and the space 
will be empty. The latter would hold if equation G 
in figure 5 was replaced by its logical inverse 
(reversing the inequality). A central characteristic 

Figure 5 

Schematic of Solution Space 



ORIGINAL PAGE IS 
OF POOR QUALITY 


(and strength) of linear inequalities filtering is 
that a given equation is either correct, or incor- 
rect (an eclipse is either correctly assigned or 
not), and if correct we will see below that the 
limits are very sharp indeed compared to the path 
width. An erroneous equation may (but need not 
necessarily) create an inconsistent sec with an 
empty solution space. 

The minimum deletion linear inequalities filter 
is defined as follows. Each equation of condition 
will have a truth probability assigned to it on the 
basis of historical soundness and observational 
characteristics in the manner of Muller & Stephenson 
(1975), the actual values used here being from 
Muller (1975). If the actual set of linear inequal- 
ities provides a solution space, it will be used 
subject to Che data consistency test below. If the 
set is Inconsistent, the group of observations with 
the minimum joint probability of truth, whose 
deletion creates a solution space, will be thro\mi 
out. This is exactly analogous to the usual 
process of deleting a few bad data in a least 
squares fit, when their residuals after the fit 
exceed a given multiple of the overall standard 
deviation of the fit (often 3*a). 

No solution by least squares or any other 
filter is any better than that which is shown to be 
sound by an adequate data cor.sistency check (redun- 
dancy). Erroneous equations will have one of three 
effects on the linear inequalities filter. [a] The 
equation is distant from the bounded solution space, 
but is oonsi stent with it (for example equation G in 
figure 5). The result: no effect on the solution or 
quoted uncertainties! [b] The equation is distant 
from the bounded solution space, but is inconsistent 
(as would be the case if we replaced G with Mi) , 

In this case, it will be discovered and deleted in 
the minimum deletion set of inconsistent equations, 
providing only that good data outnumbers the bad by 
an adequate margin (without which we have failed in 
any case). The result: no effect whatever on the 
solution! [c] The erroneous equation intersects the 
true solution space, thereby reducing its area. 

There are two possibilities. {i} The offending 
equation forms a bound near another (presumably 
correct) equation which limits the solution space. 
If, for example, equation E is wrong, there is no 
serious effect on the solution space, and the quoted 
result (centroid of the space) and error estimates 
(distances to the bounds from the centroid) are 
hardly affected. Equation E would not be identified 
as ’’wrong'* but it would have a negligible (or at 
least very limited) effect on the solution space, 

{Li} The worst possibility is illustrated by repla- 
cing equation D with ^D. Now the '’solution space” 
will be incorrectly identified as the small region 
between equations A, C, and /up. The (centroid) 
estimate of the parameters will be moved by about 
1.5 standard deviations, which is serious enough to 
avoid if p OSS ible , and the so lu tion uncertainty 
would be g ro s s ly un de re s 1 1 mat e d . This . can be 
avoided by demanding that a simple data consistency 


test be satisfied. We will simply delete any 'i>'ngle 
equation (correct or not) whose removal would 

increase the solution space sixe (as 
removing Mi surely would). If was incorrect, we 
provide a solution space very close to Uie corre.L 
one. If the correct space was the small one, we 
have provided a more conservative solution consist- 
tent with the estimated errors. Wrongly identified 
large solar eclipse observations have a tendency to 
be freest u in error (• .,/. assigning Rome insteau of 
Athens as the place of observation). Bad equations 
do not have a high probability of Intersecting the 
small solution space, class [c] above, and even if 
one does do so, application of the data convsistency 
(redundancy) test will minimize its effect on the 
quoted solution. A limited number of erroneous 
equations have a high probability of having no 
effect a'll on the solution! 

We have seen two extremely important advantages 
of this filter for the processing of historical 
eclipse data. First, the more in error an erroneous 
equation is, the lees the chance that it can have 
any effect on the solution. Second, a chance 
erroneous equation intersecting the solution space 
can have at worst only a modest influence because of 
the adopted data consistency requirement. This is 
precisely the kind of stability we need, because 
these data are either correct and very precise in 
the imposed limits, or they are very, very wrong. 

This stability is in marked contrast to least 
squares, where really bad equations of condition are 
disastrously emphasized by the least fibres rule. 

It is worth repeating that least squares equations 
of condition are only a poor approximation to the 
actual linear inequalities imposed by central 
eclipse observations, and losing sensitivity by 
using least squares is quite unnecessary (and very 
unwise in view of the relative paucity of good data). 

Muller & Stephenson (1975) undertook to eval- 
uate the equations of condition by a rather complex 
iterative scheme. This was simplified in Muller 
(1975) to a process of single parameter solutions 
using partials with respect to the observables as 
conducted below. Determination of the solutions is 
undertaken directly in AT, and since we have 
reasonable a-priori on the final solutions, it is 
convenient to express the observational bounds in 5AT: 

5AT = AT ^ '-at , ^ (24) 

observed nominal 

Partial derivatives of AT with respect to h can be 
determined by numerical calculation, thus allowing 
this parameter to be determined, while w/ulwill . 
follow from the fit of equation (18) to the bounds 
in AT versus time, and use of (17). 


Observation of Large Solar Eclipses 


A total eclipse of the sun begins quietly, a 
small nick appearing on the western limb of the sun. 
Those unaware of the impending eclipse notice 
nothing, unless clouds or a chance reflection in a 
pool of water happens to afford the right degree of 
solar filtering. The moon moves inexorably across 
the sun*s disk, taking between 60 and 90 minutes in 
the crossing. At 90% eclipse, the darkness is 
equivalent to moderate clouds and the curious casual 
observer might glance up to see what had obscured 
the sun in a cloudless sky - but he would be quite 
dazzled by the sun and would not detect the cause. 

A few minutes later, at 99% eclipse, the air would 
be chilling perceptibly, and the darkening would 
equal that of a very heavy overcast. If our 
unsuspecting observer now glances at the sun, he is 
just able to perceive the brilliant solar crescent. 
Something is decidedly wrong. It is barely 60 
seconds to sunset, at high noon! 

Day birds hurry to roost, night birds fly out. 
Animals and even insects hush without exception, 
every living creature tensely aware that scmething 
is about to happen. Suddenly, the sun is cut to a 
thin line, and shadows begin dancing on the ground 
as at the bottom of a wind-blox^ pool of sunlit 
water. All eyes turn in terror or disbelief from 
the ground, drawn to the sun. The thin line of 
crescent breaks up into incandescent pulsing beads 
of intense light, far brighter than Venus stavs; the 
beads seem to move, to corpus culate, flashing on and 
off. During the previous hour, the sun has dimmed a 
thousand times. In these last four seconds, it 
falls another thousandfold, and the umbral shadow 
races in from the west at three times the speed of 
sound. The sun WQtontly disappears completely, 
only a thin ring of pearly light visible about the 
moon*s black rim, punctuated by tiny, blood-red 
licking flames I The priest falls to his knees 
murmuring the mseri'-mei; screams of terror die 
unborn in the dry throats of thousands; the only 
sound is the flapping wings of disoriented day-birds 
flopping on the ground unable to see to fly, and 
night birds leaving for the night harvest of 
insects - 

The sky to the west is brooding black and 
darkening, while to the other compass points, the 
red, thin clouds of distant sunset reflect saffron 
light garishly into the intense faces staring 
hypnotized at the departed sun. The air chills, dew 
falls, a light breeze springs up from nowhere. The 
night darkens, and stars appear, a few at first, 
many more as the eyes adjust to the darkness. The 
sun stands now in a halo of pearly, flamelike 
streamers which seem to grOh) as the endless time 
passes - extending first to five, then ten solar 
diameters like a great celestial b lack- cent e re d , 
filamentary-pet ale d flower. The deep sunset sky now 
extends the full 360° of the horizon in all direc- 
tions. All minds are obscessed with the one thought, 
the unspeakable question: will God return the sun to 


us, or is this His final vengeance? 

All eyes are raised to the dethroned sun. 
Instantly, the night is turned to day as the first 
bead of light escapes the moon*s limb. To dark 
adapted eyes, it irradiates incandescently , throwing 
a bolt of blinding light in the pattern of a great 
cone of streamers arching down from the heavens 
directly into each bpred soul, personally. For a 
brief moment, it takes the aspect of a great ring 
set wjth a diamond of unearthly brilliance, bound 
around the black lunar disk by the thin line of 
rapidly fading inner corona; then it is too painful 
to look any more. The ground ripples again with the 
fantastic shadow shapes, and the western sky bright- 
ens rapidly, even as the shadow visibly retreats 
into the still darkened eastern horizon. Maybe - 
just maybe - it isn't the End of the World after 
all ... 

If you haven't thrown this paper in the general 
direction of the circular file, you are probably 
thinking that this description is either grossly 
exaggerated, or is a misguided attempt at literary 
prose. The latter might be true, but the former is 
not. It is difficult to convey the impressive 
spectacle afforded by a total solar eclipse with 
words; it is really impossible to do so with photo- 
graphs - though Muller (1975) references some plates 
from the literature which provide a rare hint of the 
impression. 

I have witnessed four total eclipses of the sun, 
and certainly cannot claim to have witnessed every 
wonder of nature, but there is little room for doubt 
that a total solar eclipse is the most humbling and 
impressive natural Spectacle which nature affords us 
on this planet. The twentieth century astronomer, 
and aboriginal native alike, share the same feel- 
ings - an impression of mind and soul which bridges 
the ages of man aS perhaps no other event can do. 

We are dealing in this study with a basically 
subjective phenomenon, and it is important to 
understand the subjective impression. The awesome 
impact of such an eclipse is undoubtedly the reason 
why we find so many more historical references than 
we have any right to expect. If an eclipse happens 
in a center of civilized man, it is recorded, a 
hundredfold, and this explains why we find such a 
deep historically recorded tradition in this area. 

Scientists correctly prefer quantitative proofs 
when it comes to matters of observational resolution 
and accuracy. There have been two substantial 
exp eriments which show that m trained observe rs can 
reliably distinguish the different j between a total 
and nontotal eclipse with a resolution at the path 
edge of the order 200 meters! Tlie eclipse data 
recorded and analyzed by Halley (1715) is marvel- 
ously detailed, and as reduced in Muller (1975) , 
achieves this level of resolution. Halley even 
correctly notes the conditions of the Bally beads in 
this eclipse. The franaaetiono of the Tlhiminating 
Engineering floaiety (1925.) documents an experiment 
in New York City which shox^s that a resolution of 
200 meters is achievable in practice. 
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Figure 7 


Schematic Total Eclipse Light Curve 


Typical Observational Details 



Figure 6 is a schematic of the light curve in a 
total eclipse, and is similar to the observed curve 
provided in the Transactions of the Illiminating 
Engineering Society (1925). The key to the high 
geometric resolution at the path edge is the very 
bright light provided by the smallest portion of 
uneclipsed sun, The corona (solar atmosphere) 
visible in an eclipse has a brightness of about 10“^ 
of the whole solar disk. Light scattered into the 
eclipse zone from outside it can increase this 
by about a factor of ten. In the first few km 
outside the path limit, the light level rises to 
10“^ of the sun, or three to four orders of magni- 
tude. We do not need resolution of 200m for the 
purposes of this study of ancient solar eclipses, 
and the km level of resolution guaranteed by the 
very steep lighting curve is more than adequate. 

In a typical annular eclipse, several percent of the 
sun remains uneclipsed. For this reason, very few 
annular eclipses have been adequately observed in 
antiquity. It is important to note that until the 
eclipse reaches 99%, the Solar crescent is not 
discernablc by the naked eye in a clear sky. The 
resolution of naked eye observation of a total 
solar eclipse is of the same order as modern, manual 
telescopic observations of the moon. That is why 
the ancient eclipse observations are so powerful, 
and it is crucial to use the path limits rather than 
the widths in the data reduction. 

In the reliable interpretation of historical 
records, one must be conversant with most of the 
s ub j e c tive phenoraen a in a to t ai eclipse. New t on 
(1970) and (1972) pioneered in this area, Stephenson 
(1972) and Muller (1975) carrying on. All agree 
that the statement ‘'stars seen" is not a reliable 
indication of totality. Figure 7 lists the observa- 
tional indicators adopted in Muller (1975). 

Following Stephenson (1972), Muller & Stephenson 
(1975), Muller (1975) and this paper admit an 
observation to the final solution only if it 
unambiguously indicates the precise magnitude 


Strong Indications of Totality 

T The sun disappeared 

2 The sun suddenly disappeared 

3 Darkness Suddenly fell 

4 No part of the sun remained in view 

5 Sun completely covereU 

6 Ring or bag seen around the sun 

7 The sun set in the daytime 

8 Chinese astronomer’s term chi (complete) 

Ambiguous Indicators 

9 stars seen 

10 Gi'eat eclipse 

11 Eclipse like none before it 

12 Great darkness 

13 Sign in the sun 

14 UGrrible eclipse 

Strong Indications of Nontotal ity 

15 Crescent rotated around the sun 

16 Flashing, corpus culating lights 

17 Sun appeared like a (two to five) day moon 

achieved (total or denial of totality). A statement 
of very large magnitude is acceptable for very 
ancient data (ca. -1000). 

For an eclipse record to be usable in the 
astronomical determinations, three criteria must be 
satisfied with a Sufficiently high probability of 
truth (greater than .50 in this study). They are 
listed in Figure 8. Identification (of which 
eclipse was observed) is not usually the cause of 
fatal difficulties, although there are such notable 
exceptions as the Eclipse of Hipparchus. Total 
eclipses at a given station are very rare (once per 
360 years on average), and it is usually enough to 
know the King’s relga, or a rough date, in order to 
unambiguously identify the eclipse. It xs not 
necessary to have the time of day, the day, or even 
the year, and the great strength of this data is 
that it is untimed. The question of magnitude has 
already been discussed, and the critical detail 
which is almost invariably omitted is the observer's 
location. He no doubt feels It is implied, or the 
need (to twentieth century science) for him to state 
it simply fails to occur to him! To appreciate this 
dif ficulty it is only necessary to read a few hundred 
pages f rom Newton (1970) , (1972) , or Stephenson (1972) ! 

Figure 8 

Three Observational Requirements 

I identification of which eclipse was observed 

II The precise magnitude observed 

III The place of observation (usually from context) 
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The Eclipse Observations 


Figure 9 


The historical analysis underlying the 
selection of the data adopted in this study 
conies almost entirely from Newton (1970), 

(1972), and Stephenson (1972). The sources 
include original records from the ancient near 
east including the Babylonian Astronomical Texts 
and tablets from early civilizations in this area; 
the Greek and Roman classics, upon which every 
study before Curott (1966) relied exclusively, 
but which are found to be a small source in this 
paper and Stephenson (1972); the Chinese Annals 
and Astronomical Treatises which Stephenson 
pioneered, thereby substantially increasing the 
available reliable data; and the Medieval Chron- 
icles and Texts considered by Newton (1972) and 
Stephenson (1972). Most of the observations 
noted in this paper are analyzed in Muller & 
Stephenson (1975) according to the criteria we 
adopted in our collaboration. Muller (1975) 
contains the final analysis of all observations. 

Figure 9 lists the observations admitted in 
Muller (1975) with the columns as follows: the 
date; place of observation with its latitude and 
longitude 0 to 360° east of Greenwich; observed 
phase (Total, Annular, Partial {totality denied}, 
or timed Contact); lower and upper limits of 
bounds on 5 AT as defined in (24), where the 
nominal AT is from line one of Figure 10 and 
where the coefficients A, B and C are entered 
in (18); the partials of AT with respect to n 
and U in the units of this paper; and the 
probability that the obser-vation is valid as 
assigned in Muller (1975). One of the data in 
Muller (1975) has since been ru^ed otit, thanks 
to valuable correspondence from G. J. Tooraer, 
who noted that Bagdad was not yet founded as a 
city in 693 ADI As can be seen below, data at 
this epoch plays no role in defining the mean 
solutions. It is humorous in retrospect that 
we checked the founding dates for these cities 
carefully for every data point except this one 
which was inadvertently overlooked. This is 
a p e r f e c t examp le of the law o f pe r ve r s 1 ty : th e 
worst that can happen, will happen, if you don't 
prevent it! 

Toomer (1974) p rovi des a fascinating 
argument which identifies the Eclipse of 
HipparGhus as that of -189 March 14 on the 
basis that only the elements of this eclipse 
give the lunar parallax which Hipparchus quotes as 
resulting from his eclipse reduction! Unfortunately 
it was necessary to make some assumptions about how 
Hipparchus would have s imp li fied his calcula tion , 
and some doubt can still remain about this most 
Intriguing of all eclipse identifications. Even if 
this date is taken, there is considerable uncer- 
tainty as to what part of the Hellespont 
provided the observation. It is likely that this 
data will never be sufficiently defined to be of 
use in solving for the accelerations. 


Date 


Place 


Listing of Observations 
Lat Lon 0 LL UL 




1567 04 09 

Rome 

41.90 

12.48 

T +110 

+160 

1.4 + 22 

1.00 

1560 08 21 

Coimb ra 

40,20 

351.58 

T — 

+230 

- 1 - 10 

1.00 

1361 05 05 

Mt Sumelas 

40.67 

39.64 

T — 

+ 95 

0 0 

.90 

1267 05 25 

Constant. 

41.02 

28.98 

T — 

+150 

+ 4 + 50 

.80 

1241 10 06 

Stade 

53.60 

9.48 T -160 

— 

- 2 - 30 

1.00 

1241 10 06 

Lambach 

48.10 

13.90 

T -265 

-- 

- 2 - 30 

.80 

1241 10 06 

Sheftlam 

47.92 

11.42 

T — 

+ 42 

- 2 - 30 

.80 

1239 06 03 

Toledo 

39.87 

355.97 

T + 10 

— 

- 1 - 12 

.95 

1239 06 03 

Ce rrato 

41.95 

355.48 

T — 

+276 

- 1 - 12 

.98 

1221 05 23 

Kerulen R. 

48.18 

115 . 90 

T — 

+214 

- 4 - 50 

.98 

1178 09 13 

Vigeois 

45. 38 

1. 52 

P — 

+185 

+ 2 + 24 

.98 

1133 08 02 

Salzburg 

47. 80 

13.06 

T — 

+ 70 

- 2 ^ 24 

1.00 

1133 08 02 

Vysehrad 

50.06 

14.42 

P 

+ 75 

- 2 - 24 

1.00 

1133 08 02 

Kloster.i’ 

50 . 87 

6.07 

T — 

- 82 

- 2 - 24 

.80 

1133 08 02 

Corvei 

51. 82 

9.43 

T -180 

— 

- 2 - 24 

.80 

1133 08 02 

Liege 

50.63 

5.58 

P -315 


- 2 - 24 

.70 

1124 08 11 

Novgorod 

58.50 

31.33 

T -156 

— 

+ 5 + 56 

1.00 

1079 07 01 

Alcobaca 

38. 73 

350. 87 

T — 

+ 5 

- 3 - 36 

.80 

968 12 22 

Constant. 

41.02 

28.98 

T -470 


-14 -128 

1.00 

916 06 17 

Cordoba 

37. 88 

355.23 

T — 

+280 

- 5 - 64 

.80 

522 06 10 

Nanking 

32.03 

118.78 T — 

+100 

+ 7 + 87 

.50 

516 04 18 

Nanking 

32.03 

118.78 

A — 

+130 

-13 -128 

.50 

48401 14 

Athens 

37. 98 

23.73 

T — 

+154 

+12 +157 

.50 

120 01 18 

Lo-yang 

34.70 

112.47 

P + 90 

— 

+13 +185 

.98 

71 03 20 

Chaeronia 

38.40 

22.90 

T - 62 

+ 38 

-11-110 

.50 

65 12 16 

Kuan g- ling 

32.42 

119.45 

T — 

+ 65 

-15 -150 

.80 

-135 04 15 

Babylon 

32.55 

44.42 

T — 

+125 

-16 -170 

1.00 

-180 03 04 

Ch' ang-an 

34 . 34 

108.90 

T — 

+146 

+17 +222 

. 95 

-187 07 17 

Rome 

41.90 

12.50 

T - 53 

+247 

+19+267 

. 50 

-197 08 07 

Ch' ang-an 

34.34 

108.90 

A — 

+245 

-54 -540 

.90 

-321 09 26 

Baby Ion § 

32.55 

44.42 

C -450 

-210 

-15-73 

1.00 

-600 09 20 

Ying 

30.34 

112.25 

T -100 

+760 

-26 -262 

.80 

-708 07 17 

Chii-fu 

35.53 

117.02 

T -350 

+575 

-23 -244 

.90 

-762 06 15 

Nineveh^ 

36.40 

43.13 

T +216 


+50 +650 

. 30 

-1130 09 30 

Gib eon 

31. 85 

35.20 

T -720 

+160 

+26 +365 

.50 

-1329 06 14 

An-yang 

36.07 

114.33 

X -515 

- 40 

-35-390 

.80 

-1374 05 05 

Ugarit 

35.62 

35 .78 

T +130 

1125 

100 1340 

.80 


tKlosterrath; only inconsistent observation in the set. 
^Timed contact: ertor is larger than expected, see text. 
^Was not included in solution: perhaps should have been. 

Figure 10 


# 

Run 

n 

03 /W 


A 

B 

C 

1 

Nominal 

-32.5 

-26.0 

0.000 

25. 

121. 

40.95 

2. 

3AT/3n 

-37.5 

-29. 1 

0.000 

25. 

121. 

45 .82 

3 

3AT/3Si 

-32.5 

-26.0 

-. 412 

25. 

121. 

40.95 


'cy“^ 10‘ 


r“i 


'ey 


s • cy""^ Sss • cy"* 


- 10 - 


Units 
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Two observations in figure 9 are controversial, 
and it seems that every Worker in this field manages 
to choose one or two observations which nobody else 
will touch! For -1130 September 30 we have support- 
ing analysis, Sawyer (1972) and Stephenson (1975), 
but several recent correspondents have noted 
evidence that Joshua may be dated in the 15th 
century BC. If this is the case, then Stephenson 
(1975) would prove that there is no real eclipse 
associated with Joshua and the observation Would 
have to be dropped. Some new arguments are given 
in Muller (1975) to support Plutarch's famous 
eclipse of 71 March 20. Toomer (private communi- 
cation) points out that Plutarch's reference to the 
corona occurs in the very same words in Cleomedes 
in a context drawn from earlier Greek astronomers. 
This shakes my confidence in the argument of Muller 
(1975) because the key point there is that Plutarch 
is behaving like an eye witness. It is still 
possible that he would use a particularly apt phrase 
from an earlier source, but the argument is x<7eak- 
ened. It will be seen below that these two data can 
be deleted without materially changing the solutions, 
and the remaining observation set appears very 
strong. 

The Equinox Observations 

An observation of a solar equinox (declination 
zero) given in local time for a knom place provides 
a direct measure of AT (and hence a)/w) independent 
of the moon. As shown in Newton (1970) , a small 
bias (by ancient standards) in the setting of an 
equator plane is fatal to these data unless it can 
be independently determined. As Newton notes, this 
error has the opposite effect on AT determined from 
vernal and autumnal equinox observations, allowing 
a solution for both equator bias and AT if data from 
both seasons is available from the same instrument . 
The estimate of AT from equinoxes is merely the 
difference between the observed Greenwich Time (UT) 
of the event, and the ET at which the computed 
position of the sun equals the observed position: 

ATg = ETsun - UTobs <25) 

To estimate AT in the presence of an unknown bias in 
the instrument equator, given sufficient data to 
define both vernal and autumnal values of AT: 

AT =^A(AT^-AT^) (26) 

and the equator declination bias follows from: 

AS = <AT„bs-AV-6 (27) 

Newton (1970) argues that the solstices have 
insufficient sensitivity to be useful, and they are 
not considered further in this paper. 

The equinoxes of Hipparchus are quoted in 
Ptolemy's Almagest, Fotheringham (1918) is the 
first secondary source analysis. Newton (1970) and 


Muller (1975) analyze these data, and join with 
Fotheringham (1918) in noting the peculiar problems 
associated with Hipparchus’ rounding the observa- 
tions to even quarter days. The Tropical year is 
less than the Julian year by only 0,188 hours, and 
this means that if perfect observations are rounded 
to 6 hours, there is a series of about 32 observa- 
tions at exact intervals of 6 hours, and then a 
discontinuity of 6 hours occurs as the annual 
deficiency of 0.188 hours finally reaches the 
rounding level. This curious phenomenon makes the 
interpretation of Hipparchus' equinoxes a matter of 
finding the arossomi* point in a sequence of 
observations, inferring the effective AT implied by 
the entire series. Newton (1970) notes that it is 
not correct to merely average the observations when 
the rounding tolerance is large compared with the 
observational accuracy. Muller (1975) analyzes 
these data in much the same way as Newton (1970) 
except: [1] the unpublished solar tables used by 
Newton are replaced by evaluation of the solar 
theory; [2] a more precise computation of the 
separate autumnal and vernal sequences is made; 

[3] the approximate observational accuracy is 
estimated from the observations themselves to 
replace Newton's heuristic estimate of 1 hour; 

[4] the observation made on the icpiKOS at Alexandria 
reported in the Almagest is included in the error 
analysis and solution estimate. The final solution 
is identical to Newton (1970) except for the 
slightly reduced estimated error: 

A^h " +4.20 ± 0.8 hours; epoch -\AS (28) 

To obtain w/w it is necessary to employ equations 
(17) and (18) . An approximate version of (18) would 
be sufficient, but we might as well look ahead and 
use the best available result (from the ancient 
eclipses below), equations (42) through (48), and 
us ing the last: 

w/w - -26. 5 ± 4.4 X 10”^^ yr”l; epoch -145 (29) 

The Islamic Equinoxes are analyzed by Newton 
(1970) and Muller (1975) , the latter differing by : 

[1] using an exact calculation of the solar theory 
instead of unpublished tables (which here makes a 
noticeable difference in the JED column of Figure 11 
compared with Newton); [2] assigning about one-half 
the estimated error compared with Newton. Figure 11 
indicates the datej place; longitude east of Green- 
wich , from whi eh the J UD , J uli an Uni ve rs al D at e ( UT) 
Of the observation is computed; the JED (ET) of zero 
solar declination from the solar theory; and the 
difference AT (JED - JUD) expressed in hours. Muller 
(1975) argues that observation Hi differs from noon 
by exactiy the equation of time and was probably 
calculated by an editor, that the solstice (//6) is 
likely to be weak, and that #9 is quoted as noon 
exactly but has the largest residual and is there- 
fore suspect. 

Figure 12 exhibits several defensible choices 


Figure 11 


The Islamic Equinox Observations 


# Yr Mo Dy Hr Place Lon° JUD*^ JEPt ^AT hr 


1 

830Sepl9 12.13 Damascus^ 

36.3 24476.8983 6.8961 

-0,053 

2 

830Sepl9 13.00 Bagdad 

44.4 24476.9129 6.8961 

-0.403 

3 

831Marl7 2.00 Bagdad 

44.4 24655.4654 5.5525 

+2.090 

4 

831Sepl9 19.00 Bagdad 

44.4 24842.1629 2.1351 

-0.667 

5 

832Marl6 8.00 Bagdad 

44.4 25020.7154 0.7909 

+1. 812 

6 

832Junl8 00.00 Bagdad^ 

44.4 25114.3767 4.3435 

-0.797 

7 

832Sepl8 23.30 Damascus 

36.3 25207.3646 7.3878 

+0.557 

8 

844Sepl8 21.40 Bagdad 

44.4 29590.2629 0.2882 

+0.607 

9 

851Sepl9 12.00 Nisabour§ 

58. 8 32146.8313 6.9797 

+3.554 

10 

882Sepl9 1.15 ar-Raqqah 

39.0 43469.4342 9.4911 

+1. 366 


tjUD ~ 2 ,000 ,000 days; JED is the same day as JUD. 
^Questionable data, see text. 


Figure 12 

Islamic Equinox Means 

# Solution AT hr A6 hr Comment 

1 Mean of all: = wts +.806 Newton’s Preference 

2 All: //6 wt = ^ 'f.937 Newton’s Preference 

3 Delete: 1,6,9 +.766 Muller's Preference 

4 Bagdad set +.708-1.243 Removes Equator Bias 

5 //8 Alone +.607 Instrument Reset? 

6 Mean of 7,8,10 +.843 Late Observations 

Adopted Mean: +.778 ±.140 

for averaging subsets of the data to find AT, 
including the choices of Newton (1970) and Muller 
(1975). Newton also comments that a value near 
0,7 hours could be taken by any reader as reasonable. 
The adopted mean and standard deviation adequately 
encompass the alternate data choices, and lie close 
to the preference of Muller (1975) on line 3. The 
better data are given to 1 hour, which would imply 
a standard deviation from roundoff of about 0.33 hr. 
It is likely that 4 or 5 of the observations reach 
this level of observational accuracy, and the 
adopted result is: 

ATj =+ 0.778 ±0.14 hours; epoch 840 (30) 

Converting to 0)/ai as with (29) yields: 

OJ/OJ^ - 22.2 ± 2.8 X 10"^^yr“^; &pooh 840 (31) 

Newton (1970) provides an analysis of the 
stated mean positions of sun and moon given by the 
so-called Hakemite Tables. Muller (1975) also 
analyzes this presumably ephemeris position: 

Xg - 254. 7658° ; X^^ - 270. 6868° ; X^^ - Xg - 15 .,9210° (32) 

Epoch: lUo on ^ Cairo Me an Time 30 November 1 0 0 0 AD . 


Newton (1970) rather arbitrarily, in the 
opinion of Muller (1975), assigned an uncertainty of 
30 arose con ds to this ’’observation”. This uncer- 
tainty might not be too optimistic for the obser- 
vational mean; but mean with respect to what? The 
steps necessary to obtain a true mean ephemeris 
position on the abstract equator are subtle even for 
modem astronomy students. 

Did the framers of this ephemeris first make 
a long series of observations and then set the 
instrument equator before making further ephemeris 
source observations? Did they write their ephemeris 
with respect to the best observed value of the true 
equator, or were they satisfied with an instrwvent 
ephemeris? Did they recognize the difference? We 
have seen evidence from the Islamic Equinoxes, 
figure 11, that the Bagdad equator was reset with 
improved accuracy. It is therefore probable that 
the Cairo instrument would also be reset at various 
times before, during, or after the observations were 
made which yielded the stated ephemeris position. 

We could perhaps rely on their equator if we knew 
that they solved for the abstract, ephemeris equator 
from the data, but I am unwilling to assume that they 
did this. It seems generous to assume that the Cairo 
instrument and/or ephemeris equator was three times 
as accurate as the computed equator bias at Bagdad 
near epoch 832 (line 4 of figure 12). This yi&lds 
the estimate of error for the AT implied by the 
Hakemite Tables as computed by Muller (1975): 

AT.J, = 3223 ± 1500 s; epoch 1000 (33) 

and converting to w/o) as in (29) yields: 

a)/03^ = - 32.2 ± 13 X lO'^ yr”^ ; epocTz 1000 (34) 

The weighted mean of equations (29), (31) and 
(33) provides the best estimate of w/w independent 
of the moon from the Timed Solar Data: 

Ui/Ui = - 23.8 ± 2.3 X 10”**yr"^ (35) 

TSD 

Since, as noted above, there is a very tight 
linear relationship between the parameters ii and 
(L/w, we can compute the former value corresponding 
to the latter in (35). Again, looking ahead for the 
most precise equation (47) we find: 

^TSD ~ ~ 2.9.1 ± 3.7 "cy~- (from 

Solutions from the Large Solar Eclipses 

The observations accepted for inclusion in the 
statistical regressions have been tabulated in 
figure 9. Determination of the solutions for n and 
0 )/ w will be undertaken in the single parameter AT by 
using the partials given in figure 9 by way of 6AT 
as defined in (24) whei*e the nominals come from line 
1 of figure 10. The 6AT upper and lower limits for 
any other parameter or coefficient choices can be 
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found from: 


6AT 


nom 


.. 3AT ^ 
6?2 + 
9n 


6U 


9AT 


5ft 9AT 


- ( 6A + 5B*T + 5C*T^) (37) 

where the parameter changes are in the sense new 
value minus nominal (6n = auxi (37) can 

be used for both upper and lower limits as needed. 


Figure 13 plots the 6 AT upper and lower limits 
for the observations of figure 9, mapped to h--35 
and -31 ”cy"^ via (37). The limits corresponding to 
these two accelerations are marked on those bounds 
which differ sufficiently (see the Key also). From 
this plot alone, or the equivalent calculation from 
figure 9 and (37), the determination of n and a 
complete AT polynomial can be made. The plot is 6AT 
versus epoch of observation. 


Figure 13 


Observational Bounds 



H 


r> (*i cn fo 

I >— H 1 



1567 1130 

70 72 


^ -H 1 ^ 

c>i o 

Max /5AT 


AT range with limit* 

AT limits at n = -31, ^-35 "cy~^. 

instantaneous (and Modern) AT, 

Epoch (astronomical year) of the observation. 

At values from Newton (1970,1972) at n - -35. 

Final fit to mean AT polynomial. 

Timed Contact , increasing delay : n- -35; minutes . 

Maximum departur e of AT ins tan taneous from mean 
between two epochs (read difference from abscissa 
scale above curve). 



Figure 14 


The solution space in AT expressed as a poly- 
nomial (18) is the set of all polynomials which can 
successfully pass between the limits shown in figure 
13. Inconsistent observations (if any) are deleted, 
and medium/short term AT variations (few centuries 
time base) are smoothed. The solution space in fi 
is defined as the range of n over which there exist 
polynomials for AT consistent with the data and the 
minimum deletion data set required by the filter. 

The solution for A is undertaken first. 

Considering the curve ’’instantaneous AT” in 
figure 13 indicates that the only identifiable 
inconsistent equation is that from Klosterrath, 

1133. Such a recent observation has no significant 
effect on the mean solution. Striking a rough mean 
of the modem data (since 1000) in figure 13, and 
turning attention to the ancient observations which 
dominate the overall determination of the acceler- 
ations, it is clear that there exist AT polynomials 
for the bounds computed at A = -35”cy"’^. It is also 
obvious that for w differing sufficiently from this 
value, one or more observation bounds will make it 
impossible to fit the AT polynomial between the 
limits. Each violated limit has a probability 
associated with it, and as we try to force n to 
wider values, thereby violating further limits, 
the probability that the given value of h is ruled 
out grows according to the ordinary laws of 
probability. To demonstrate the bounding of n it is 
only necessary to find the first few paij^s of 
observations which stand as upper and lower bounds 
to the parameter. 

The dashed curve in figure 13 labeled AT” 

satisfies the observations at n - -35”cy”^ , and this 
is the place to begin in constructing the solution 
bounding pairs shown in Figure 14. At -31”cy"‘^ it 
is impossible to fit AT, and several observations 
operate against this value (for example 65 and 120). 
Computation reveals that the cutoff point for this 
pair is h- -31. 9 "cy“^ , and this value is entered 
as a bound on line 3 of figure 14. 

Figure 14 is cons tructed by f inding the 
s e quen ce of upp e t an d 1 owe r b oun ding p ai r s as just 
discussed. The columns include several direct and 
us ef ul numbers . [1 ] Th e ob se r vation p air numb e r . . 

[2] The observation pair epochs. [3] The differ- 
ence between the appropriate UL and LL of figure 9 
which expresses the ’’width” of the common area 
through which the AT polynomial can be passed, for 
n - -31”cy“^ . Figure 9 is computed for -32 . 5”cy“^ 
as noted in figure 10, but thn node acceleration 
correction at this w from (14) and (9) gives 
a correction to the apparent of about +1.5 ”cy~^, 
and th is is added in to the solution a t this point . 

[4] This is the difference of the partials 9AT/8« 
for the two observations in the pair from figure 9. 

[5] The bound is computed, therefore, from: 

^bd " / 3AT. Exact computation, and the 

need to account for the change in the nodal acceler- 
ation correction, makes the numbers given in figure 
14 differ slightly from rounded values in figure 9. 


Solution Bounding Pairs 


# Observations 6t 3AT n Estimate 


Upper Bounds 




1 -135, 120 

-35 30.2 

-29.8 .98 .99 

-36.7 t 2.5 

2 -1374,-1130 

-30 77.5 

-30.6 .40 .97 


3 65 , 120 

25 29.3 

-31.9 .78 .96 

-37.4 ± 2,7 

4 -1374,-1329 

170 141. 

-32.2 .64 .82 

-3i.3 ± 1.9 

5 71, 120 

52 25.0 

-33.1 .49 .49i 

-34.4 ± 2.0 

Lower Bounds 




6-1130, -600 

260 54.4 

-35.8 .40 .40i 

-34.4 ±2.0 

7 -600, -180 

246 45. 0| 

-36.5 .76 .78 i 

-34.3 ± 1.9 

8 -180, 71 

208 29.2 

-38.7 .47 .86 ' 


9 -708, -180 

496 41.8 

-42.9 .85 .95 

-37.4 ± 2.7 

10 -1374,-1329| 

1770 141. 

-43.4 .64 .98 


n -1374, -708| 

1605 129. 

-43.6 . 72 .99 

-36.7 ± 2.5 


The probability associated with a given pair- 
bound is clearly the product of the truth probabil- 
ities for the two observations , since the inference 
fails if either is wrong: 

= P(t)i*P(t)3 (38) 

In the example, (.80 X .98 - .78) as indicated in 
the column P^^ of figure 14. Other bounds are listed 
in mono tonic order in figure 14. 

Moving from the middle of figure 14, up and 
down through the upper and lower bounds, the 
probability that w itself has been bounded 
increases. As each successive bound is crossed, the 
individual probabilities for each bounding pair P^ 
will compound in yielding the overall inference 
probability P^^^ that the actual value of ?V is 
within the solution space bounded by that value. 

If each bound arises from independent observations, 
th en the p rob ab i 1 i t y th a t h is bo un ded ins i de t he 
nth P^ will be ; 

^inf = Lvfl (39) 

- : i= i 

Reading up oi* dowm column P^j^£ in figure 14 from the 
middle will provide the running product (39). It 
will sometimes happen that a pair of bounding 
Obse rvations includes an obse rvat i on whi ch has 
already b e en in el u de d in a p re vi o us pair . It would 
be incorrect to include the same observation's 
probability of truth twice in (39)i Consider two 
bounding pairs with three observations p, q, and r: 



LSE 
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p with q, and p with r. This pair of crossed bounds 
must be replaced with the single bound: p witl) 

(q or r) in evaluating (39), and this is reflected 
in figure 14. This process is much easier to do, 
than it is to explain. 

The final column of figure 14 demonstrates the 
conversion, at last, of the probabilities into the 
more familiar solution with uncertainty. This can 
now be done without the devastating loss of sensi- 
tivity which would have accompanied the expression 
of equations of condition in least squares form. 

For any pairs of lines in figure 14 where 
values from upper and lower bounds are roughly 
equal, we can form a weighted mean solution for 


The estimate of w/w comes from (42) and (17): 


= - 24.3 ± 2.0 X 10“'*yr~ 
ad)/03 » 0 . 2 i-f i i?. ' >1 fi ) 


( 43 ) 


The various linear combinations of the para- 
meters can be formed to provide very useful conver- 
sion relationships. Since the A and B coefficients 
in (42) are relatively constant over n, it is 
sufficiently accurate to express the linear relation- 
slUp to C from (17) and (23): 


is. 

3n 


.980 


(H) 


h (Estimate) from the corresponding to the 

lines. An approximate standard deviation of this 
mean can be computed from Pj^^f difference 

between the values, to make the probability 
roughly V3 that the true value lies within the 
standard error. It is unnecessary to make this 
calculation with great precision, and an exact 
formula would be rather complex, and in any case, 
require assumptions about the Gaussian distribution. 

Muller (1975) prefers the n estimates near P^^^ 
probabilities between 1 and 2a (i.t?. .67 and .95), 
rather than the remote limits at the higher 
probabilities; nearer bounds to the crossover point 
in the middle of the figure are probably most 
significant. An n = -34.5 ± 3 "cy“^ was adopted, 
but values between -34.3 and -37.4 could be defended 
on various grounds. Deletion of the observations at 
71 and -1130 would tend to move n nearer -35.5, 
indicating the small effect which, arises in deleting 
these possibly controversial observations. There are 
other bounding pairs ready to enter the figure 14 
at the top and bottom to replace deleted data from 
the table. It is unlikely that any defensible 
choice could move h more than the standard error. 
Applying the correction to h from (12) , and 
retaining the Muller (1975) estimate from figure 14 


of -34.5 


^LSE 


cy 


30.0 ± 3.0 ’'ey”- 


(40) 


To have a consistent and accurate ephemeris, it is 
necessary to apply the above noted corrections to 
the ILE (1954) as follows: 


The result (42) can be mapped to other values of 
if desired. From (40), (42) and (44): 


C = 


.980^ + 8.90 


(45) 


where (as elsewhere) the units are as stated above. 
For w/w given we use (45) and (17) : 


dj/w = + .622n - 5.65 + 0.2 
and h given w/o) similarly: 

+ 1.61w/u) + 9.10 ± 0. 3 


(46) 


(47) 


Also from (17) and (42) we can obtain oi/oj given AT 
at an epoch T: 


r 20-AT , xxH 
oj/w = .635 { " 


114 

T 


aAT >> 
t 2 J 


(48) 


we have from the large solar eclipses: 


Constancy of the Accelerations 

Newton (1970), page 644, concluded that the 
acceleration(s) changed substantially near epoch 
700 and were far from constant during the historical 
period. He defined a parameter D" to measure the 
constancy of the linear relationship between h and 
m/u as follows: 

D'’ = n -* 1.6073 0)/w (49) 

which is equivalent to equation (47) above, D" 
corresponding to the constant term, and we have: 


^2 


+4.39 ”cy" 


and Q - - 0.120 


cy 


(41) 


+ 9.10 ± 0.3 ’’ey' 


-2 


(50) 


Equations (40) and (41) if applied to the lunar 
theory will provide the best fit to the ancient' 
solar eclipses. 

The best fit to AT is indicated by the dashed 
curve in figure 13, and corresponds to: 

AT = 2D + 114*T + 38. 30-T^ ±0AT=.3T^+5O (42) 

@ nt = -30.0 "cy-=' 

The uncertainty was estimated in Muller (1975) on 
the basis of deleting data and bending the AT 
polynomials, and applies to AT ^tiven n. 


Newton found various values between ±20 from 
selected subsets of his data and concluded that 
the accelerations had varied. Muller (1975) shows 
that the typical short and medium term variatioris 
in AT such as those shown in figure 13 and as 
discussed by Morrison & Ward (1975) require the use 
of at least 800 year time bases in the estimation 
of whether the parameters have varied durihg 
historical time. Periods of 400 - 800 years are 
usab le if errors larger than the parameter un cer- 
tainties (40), (43) are acceptable; less than 400 
years will be quite unreliable in this connection. 


Figure 15 


Various Determinations of D" 



Source 


Epoch Range 


1 

Newton 

1970 

-500 

+500 

+3.0t 

2 

Newton 

1970 

+500 

1300 

-6.lt 

3 

Newton 

1970 corrected 

-500 

+500 

+8.6 

4 

Newton 

1970 corrected 

+500 

1300 

+7.4 

5 

Newton 

1970 p. 641 

1650 

197a 

-18.0^ 

6 

Muller 

1975 

1567 

1975 

+10,5 

7 

Muller 

1975 

various 

+5 to +12 


Muller 

1975 & This paper mean 

+9.1 


At modeling error - see text. 

^ Insufficient data base. See note page 23 below. 

Figure 15 summarizes the debate regarding 
constancy of the accelerations. Newton (1970) 
modeled AT according to (16) in computing w/w via 
(17). to examine what Newton (1970) would have 
gotten x<rith the correct expression (18), it is neces- 
sary only to convert his (33/ w results back to a value 
of AT via (16), and plot them directly in AT. This 
was done in figure 13 near the years +1000 and -200 
and marked "70” for Newton (1970). These AT values 
can now be converted correctly Into oj/o) via (48), 
which when combined with the corresponding value of 
n, yields T>" from (49). The first two lines in 
figure 15 give Newton* s P" solutions, the next two 
lines give the solutions corrected for modeling AT 
as above. Line 5 was his speculation on the modem 
value from occult at ion data. Muller (1975) examined 
the modem occultation data of Martin (1972) and 
Brouwer (1952), including the absolutely reliable 
eclipse observation by Clavius in 1567, and found 
the D" value on line 6; other subsets of the data 
(including solid and dotted curve in figure 13) 
provided a range of solutions as on line 7. The 
consistency between the correctly modeled determin- 
ations of D" in figure 15 proves that there has 
been no significant change in the relationship 
between the two accelerations during the historical 
period. 

The question can be asked, would we fail to see 
a propox*tion(xl vartatlon in the accelerations when 
examining the constancy of D"? this test will 
indeed be weak if the accelerations vary together 
according to their mean partials (49) , If the tidal 
couple were to ch an ge , as has been s ugges ted by 
Newton (1970), then the parameters are related by 
the functions (19) and (20) , which differ from (23) 
and a change in the tidal couple would be visible 
with substantial sensitivity. The medium term 
variations in AT, figure 13, bear further study. 


I have just received a copy of Newton (1976), 
and therefore add this note in proof. Newton now 
accepts solving for constant and linear terms in 
(18), although he says that this must not be done by 
correcting the ephemeris. He fails to see that (18) 
and (3) are the same except for a constant scale 
factor, and it does not matter which is corrected. 
Although he still believes that the accelerations 
have changed, it is clear that the approach taken in 
Muller (1975) and above is now accepted by Newton. 

To quote from his page 12; "Similarly, when we deal 
with the acceleration of the sun, we should not 
alter the constant and linear terms in the ephemer- 
ides in response to an acceleration estimated from 
ancient data, since we have strong reasons to 
believe that the accelerations have not been 
constant since ancient times. Instead, we should 
estimate the constant and linear terms entirely from 
recent data, We should also estimate the acceler- 
ations from recent data and compare with the values 
estimated from ancient data." The approach taken in 
this paper xa pvecisety what Newton suggests. It is 
clear from figures 13 and 14 that the ancient data 
alone determine the accelerations, the medieval-to- 
modem data span provides coefficient B, and the 
current epoch establishes A in (18). The solar 
system integration, Muller at. (1976), found 
that the constant and linear terms in the hmax* 
ephemeris (1) do not require significant corrections. 
The entire scheme is therefore consistent and will 
correctly estimate the accelerations. 

Measuring the constancy of the accelerations is 
a basic goal of the recent studies in this field, 
and it is i mp ortant to be quite specific in c le a r ing 
up the misconceptions. Newton (1976), pages 11 & 12, 
appears to concentrate on three points of evidence 
in support of hiS contention that the accelerations 
have changed during historical time. [1] Apparent 
changes in D" (or the equivalent) are indicated by 
ancient observations near epochs +700 and +1300. 

[2] The coefficients A and B in (18) must be solved 
from recent data, by which he sometimes seems to 
mean "not ancient" and at other times ’’very recent" 
i.e. the last 2O0 years. [3] He suggests comparing 
accelerations determined from ancient and modern 
(1700-1975) data. Each of these three points will 
be considered in turn. Newton's earlier works used 
D" j equation (49), to test the constancy of the 
accelerations. Newton (1976) refers to , the 
lunar acceleration on solar time, but this differs 
from D" by only the small value of the solar accel- 
eration Vg V which is near -2,5 "cy~^: - D" + Vg' . 

[1] Having given the quoted statement above, 
Newton (1976) still employs (16) instead of (18) in 
estimating u)/w, which in turn is central to D" and 
V ' in th e t e s t f Q r c on s t an cy . Ne w t on ' s con cl us ion 
that the accelerations changed near +700 and +1300 
is based on the D" differences between the data sets 
noted in figure 15, lines 1 & 2. When his results 
are corrected, the ancient and medieval data show 
little variation, lines 3 & 4. This proves that it 
is Newton’s reliance on equation (16) Instead of (18) 
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in the presence of significant error in the rate of 
Newcomb’s solar ephemeris (3) which gives rise to 
the apparent, but quite opuPtouB change in the 
accelerations near these epochs which he notes. 
Another Way to see the problem is to note that if 
(16) is to be used, we must have T = 0 be the epoch 
at which coefficient B in (10) is zero. As Muller 
(1975) notes, it can be seen that this occurs near 
the year 1800, not 1900 as assumed (in effect) by 
using equation (16) in Newton’s works. See (62) and 
item [3] below, for examples. The solar ephemeris 
condition of equation (3) could be mcq^ped to any 
year, say between 1770 and 1800, since the epoch is 
quite arbitrary. Had Newton used that epoch, he 
would have found the 0" values given in figure 15 
lines 3 and 4 instead of lines 1 and 2. Unfortu- 
nately, Newton (1976) notes that it is all right to 
solve for A and B in (18), but fails to apply the 
implications of what is in fact the necessity of 
doing so to his current and earlier results. He 
therefore persists, erroneously, in his belief that 
his earlier data analysis indicates a real change in 
the accelerations near +700 and +1300. 

[2] As noted in the first paragraph of this 
note, the procedure of figures 13 and 14 determines 
coefficients C, B and A of equation (18) with 
ancient, medieval, and modem data respectively , and 
independently. If a determination of these 
coefficients Is desired from a more recent data span, 
the recourse must be to the next paragraph. 

[3] Newton (1976), page 12, says that Muller & 
Stephenson (1975) ignores the difference between 0" 
lines 5 and 8 of figure 15 (in effect), thereby 
missing the change in the accelerations. Muller 
(1975) investigated this last vestige of Newton’s 
argument by examining AT from 1567 - 1975 as follows: 
{a} Clavius' observation of a total solar eclipse in 
1567, AX = +104 ± 25s; {b} Martin (1972) mean at epoch 
1640, AT-+22± 5; {c} from 1650 - 1900. AT=^ 0; {d} AT 
then rising smoothly to +45s near 1975. Almost any 
D" can be obtained by fitting various data 1650-1975 
as done by Newton. A representative 1567 - 1975 fit 
results from solving (18) from three points : 1567 
(+104), 1654 (+2), and 1970 (+40). This provides 

AT = -24 + 68T+ 32T^ ^ and hence a)/u= -20. 4 X 
from (17) cor responding to the ephemeris A = -22 . 44 
’'cy“^ used in the computation. From (49) we find 
D" -+10.5 ”cy“^ as shown on line 6 of figure 15. 
Other data selections in this interval provide O" in 
adequate agreement with figure 15 lines 3, 4, and 8; 
not line 5. This is in better agreement, in fact, 
than we have any right to expect given the relation- 
ship between 6T (the time base of the data in cy) 
and the uncertainty in D" arising from medium term 
AT variation, Muller (1975): 2cy, 19; 3cy, 13; 4cy, 
8. 8 ; 8cy , 3.9; 10 cy , 2.4 ’Vcy“ ^ . Newton overlooks 
this , and his final point vanishes . 

For ait data spans which ore long enough, D" 
remains in very good agreement, observation ally 
demons trat ing th at th e ac ce le r a t i on s h ave remaine d 
constant within uncertainties (about ±3 "cy”^) 
corresponding to averages over 8- 10 centuries. 


The Earth*s True Nontidal Acceleration 

Under the assumption that G/G-0, we can solve 
for the earth’s apparent ^nontidal rotational 
acceleration from (40) and (43) via (22) : 

'^/^ANT ” IT* 3 i 4 X 10“-^yr“^ (51) 

This is about twice the currently suggested expec- 
tation from postglacial uplift, Kaula (private 
communication); see Dicke (1969) and O’Connell 
(1971). Obviously, if real rather than apparent, 
this acceleration must arise from changes in earth’s 
moment of inertia or interchange of angular momentum 
within the earth, with a time-constant long compared 
with the historical period. Yukutake (1972) 
suggests that maintenance of the changing geomag- 
netic dipole moment might account for +5 X 10“ yr“V. 
The current westward drift of the earth’s magnetic 
field might be revealing the results of an acceler- 
ating couple inside the earth, but the required 
interchange of momentum demanded by 3000 years of 
this seems rather large. Lyttleton (1965) considers 
the possibly controversial idea that the earth’s 
core is a phase change in mantle material, and makes 
some predictions of planetary conditions in the 
solar system (1969) which have proven interesting 
in view of recent observations by interplanetary 
spacecraft. With this hypothesis he resolves the 
Jeffrey’s tidal inconsistency, and finds a 
consequent nontidal acceleration of the earth’s 
rotation, Lyttleton (1976), of: 

= + 8 X 10”“yr"^ (52) 

This is in adequate agreement with (51) , constitu- 
ting another interesting prediction of the core 
phase change model. It appears that the geophysical 
theories are in a state of flux, no one of them yet 
proven. A significant &/G remains a possibility , 
however, see below. 

Lambeck (1975) undertakes a detailed study of 
the theoretical earth-moon tides , and obtains 
reasonably good agreement with observations from 
his computation, finding: n - -35 ± 4 ”cy“^, figure 1. 

Cosmology and Final Equations of Condition 

The earliest suggestion that G/G could be 
observed in solar system astronomy is due to Dicke 
(1957). If G has a rate of change, then when 
observed on an invariant time scale such as Atomic 
Time (or the rotating earth) , all planetary orbits 
expand and the periods lengthen proportionally: 

2*&/G = Ag / n (53) 

where is the orbital longitude acceleration 
ob se rve d on in Var i an t time . In the Di r ac (19 73) 
cosmology, it is the atomic unit of length which 
is changing, but the observational principle remains 
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the same except for a sign ambiguity and scale 
factor change (the 2 in equation (53) becomes 1). 

Dicke (1966) noted that if the earth’s rota- 
tion can be calibrated for all geophysically 
induced torques^ then it can be used as the cosmo- 
logically invariant time. Since the earth’s 
observed rotational acceleration cb/o) is measured 
with respect to the sun’s position, it will appear 
to be accelerated in the presence of G/G: 

6w/a>obs “ 2.0 G/G (54) 

If the earth was a completely rigid body, only 
equations (54) and (22) combined with suitable 
observations, would be required to separate 
(the true nontidal earth acceleration) from G/G, 
thereby ob€eri>aHonally determining both. This was 
undertaken in Muller (1975) with encouraging 
results. Two problems arose in this connection. 

The first, noted by Muller (1975), was to eliminate 
the necessity of solving for the lunar nodal rate 
from these observations. This problem is now 
resolved. The second problem arose in later 
analysis. 

Murphy & Dicke (1964) estimated that as G 
decreases, the earth expands due to the relieved 
compressional stress provided by gravity according 
to the formula: 

6w/t0o^g = + 0.20 G/G (55) 

This equation was reviewed with the same result by 
Nordtvedt & Will (1972). The bulk modulus of 
interior portions of the earth Is known accurately 
from seismic wave observations, and the references 
found (55) on that basis. Lyttleton (private 
communication) has pointed out that the reference 
assumes the core is incompressible, uses a model for 
the mantle which is too stiff, and ignores the known 
phase change at the intra-mantle boundary. \*Jhether 
one accepts Lyttleton 's (1965) model of the c ore- 
mantle boundary as a phase change or not, it appears 
that the expansion ratio found by Nordtvedt & Will 
(1972) is too small. An exact solution to the 
problem appears to require some effort , but approx- 
imations indicate increases of between 2*5 and 5 in 
the coefficient of G/G given in (55). Pending a 
review of this expansion question, a parametric 
study of this will be substituted for (55). Define 
e as the expansion coefficient and write: 

doj/to^^g = e-G/G ; 0.2 < £ <1.0 (56) 

where £ probably lies in the range shown. 

To include Che effects of 0/G in the recon- 
ciliation of all known effects in the earth's 
rotation, where is the true geophysical 

acceleration, we modify (22): 

^ 1.1 - (2.0 - e)G/G + (57) 


Dicke (1966) calibrated the earth’s rotational 
acceleration by way of a detailed theoretical 
consideration of the known geophysical mechanisms, 
thereby finding It is proposed here to 

solve for ft/G, and directly from the 

observations. The ancient eclipses even with the 
timed solar data above are not enough by themselves 
Co successfully solve for the three parameters. 

The full observational equations of condition from 
all applicable astronomical data will be written, 
including the optional constraint 5/0 « 0. 

Figure 16 


Equations of Condition 



Figure 16 displays the adopted equations of 
condition from various observation Sets. The first 
block gives the equations of condition which hold in 
the presence of the so-called cosmologies, 

in c 1 ud in g Hoy le & Na r lika r (1974), B r ans & D i ck e 
(1961), Peebles 5 Dicke (1962) and others. The 
second block applies to the Dirac (1973) cosmology. 
The sign on the coefficient of G/G has been set in 
agreement with his multiplicative matter creation 
postulate, and a negative value solved for G/G 
will support this postulate, whereas a positive 
value would point to additive creation. 

Line one of figure 16 comes from (40), with the 
Small correction to apparent w arising from the 
effect of G/G on si from (9) and (14). Two is (43) h 
(57), where w/oj is also subject to through its tie 
to h as in (46) . Line three is the timed solar data 
(35) , also using (57). Line four is the result of 
Morrison & Ward (1975) who reworked the Hpencer 
Jones transits of Mercury (modern data 1670-1975) . 

As already noted, timed lunar data is insensitive 
to error in nodal rate. 
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VanFlandem (1975) and (1976) approaches the 
problem directly by making observations on Atomic 
Time, where from (53) we can write: 


Figure 17 
The Solutions 


/n - /Zt /w 2'G/O (58) 

Ills most recent value (1976) for the lunar occulta- 
txons 1955-first quarter 1976 is: 

= - 36.0 + 5.0 ”cy^^ (59) 

which gives rise to line five of figure 16. Line 
six merely provides the constraint &/G=0, which will 
he used as an alternative cosmological assumption. 

In the Dirac cosmology we must replace 
equations I, 2, 3 and 5 as shown in the bottom of 
figure 16. Since it is the atomic unit of length 
which is changing in this cosmology, there is no 

dynamical expansion of the planetary orbits, 
and there is no effect on the node or on earth 
expansion. 

Planetary and lunar ranging data will probably 
yield estimates of parameters with usable sensitiv- 
ity in determining (;/0 within a few years, but they 
are not quite there yet - see for example Reasenberg 
& Shapiro (1976). The paleontological data of 
Pannella (1972) can yield estimates of d)/03 and 
Kaula & Harris (1975). It is my view that: [i] the 
uncertainties in the recent udta (ea. -50 X 10^ yr) 
are too large (they quote n ~ ~5P ± 15 *’cy”^) ; the 
data sample is limited and a very small systematic 
error in the fossil counts will obviate this result; 
[ii] the tidal couple may have changed significantly 
in even 50my, and almost certainly has over 450my 
(at which epoch they quote fi = -18±5, error estimate 
mine). Newton (1968), figure 1, finds from near- 
earth satellites, but I am concerned about the short 
data span and complexities in the theoretical 
modeling, which (if incomplete) will seriously 
impair the results. These potential observational 
results are not included in the regression, for the 
reasons stated. 

Solutions from all Astronomical Data 

Figure 17 provides the solutions from the 
equations of condition (figure 16) shown in column 
two. Line one is the so-called primitive cosmology > 
line two on Dirac's cosmology , line three assuming 
fc/C-0. Lines four and five parameterize e, (56) & 
(57) , for the primitive cosmology. Line six solves 
the Dirac equations with the Hubble rate from extra- 
galactic data (68) added, to get the final 

The value of remains between -27.2 and -28. 1 
for all cosmological assumptions. The value from 
modern observations due to Morrison and Ward (1975) 
in figure 16 line 4 is probably less subject to a 
systematic error than the ancient eclipses (line 1) 
because the latter does depend on a smaller data 
sample than one might prefer. The solution 
consistency and weights on each data set seem 
reasonable, however, and the overall solution is 


n 

Eqns 


(Vo 0 

W/W'I'NT 

e 

1 

1-5 

-27.4 1.6 

-2.3 1.5 

+4.0 4.1 

0.2 

2 

la-5 a 

-27.2 1.7 

-5.1 3.0 

+3.1 4.4 

any 

3 

1-6 

-28.1 1.6 

O.G 

+9.2 2.5 

any 

4 

1 - 5 

-27.4 1.6 

-2. 3 1* 5 

+4.7 3.8 

0..5 

5 

1 - 5 

-27..4 1.6 

-2.3 1.5 

+5.9 3.2 

1.0 


la* ia*^ 

-27.2 1.6 

-5.6 0.7 

+2.5 2.5 

any 


"cy"^ 10~^^yr^^ 10^^^ 

twith eqn. (68) added. 


probably the best choice for a final value: 

- - 27.2 ±1.7 “cy--^ (60) 

The agreement between these determinations of 
and the values of (b/o) from figure 16 line 2, 
dependent on the lunar acceleration, and line 3, 
independent of the moon, builds confidence in the 
elimination of the Spencer Jones anomaly, and in 
the constancy of the accelerations during the 
historical period (within uncertainties) as argued 
on other grounds above. 

Only the linear combination of parameters found 
from the ancient data can effectively provide the 
apparent value of w/o eorre^f^-:md:-ng to this and 
from (60) via (46) : 

fl)/u = - 22.6 ± 1.1 {0.2} X 10~“yr"^ (61) 

The parenthetical uncertainty applies if is 
perfectly known. This result remains in very good 
agreement with the value from timed solar data 
independunt of the moon (35). 

To be of practical use in the computation of 
astronomical ephemerides during the historical 
period, equation (60) must be accompanied by the 
complete expression for AT, replacing the values 
on page 87 in the Explanatory Supplement (1961), 
Using (60) , (42) and (44) we find: 

AT = 20 + 114'T + 35.55-r- ± 0AX== .3T- +50 (62) 

To have a consistent, and correct, lunar ephemer is, 
it is necessary to add two corrections to the Imar 
node expression coefficients in (2). First the 
correction to the nodal rate detertnined from a 
direct integration of the solar system: 

ft = + 4. 39 "cy"‘ (63) 

Second, the theoretical eotrection to nodal acceler- 
a t i on res ul ting f rom the tid ally in du ce d e xp an s i on 
of the lunar orbit from (60) & (14) ; 
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(64) 


Rate of Change of G 


where h this value Is added to the coefficient of T^* 

The ephemerls based on equations (60), (62), 

(63) and (64) will provide the best fit to all 
astronomical observations, ancient and modem. This 
will apply equally well to all observations made on 
Universal Time regardless of the cosmological 
assumption regarding G/G since (62) is the oboct*VCd 
apparent AT, and net S corrections arising from 6/G 
are very small. 

The values of fall in the range 4*2.5 

±2,5 to 4*9.2 ± 2,5 for the various cosmological 
assumptions. These values are in the range which 
current theoretical analyses view as reasonable, and 
it appears that there is no serious observational 
versus theoretical discrepancy - a significant 
improvement on results such as Newton (1970) who 
obtained 4*23 X 10“’^^yr“^. The issue now appears to 
be clearly drawn. If there is a real or apparent 
G/G, then there is no significant To 

geophysicists committed to one or more of the 
several theoretical explanations for a significant 
nontidal earth acceleration, this may possibly 
cause them to believe that G/G-D! Conversely, 
most CQsmologists will probably argue that every- 
thing is tidy and consistent - a reasonable G/G and 
an earth* s tidal rotation satisfying the conservation^ 
of angular momentum! There is always the poswibil- 
ity of a bit of both, as Urey frequently and wisely 
points cut. In this case, we would need a cosmology 
demanding a small fraction of the Hubble Constant 
as its rate (see below), and I am not aware of one 
in the current literature. It will be apparent 
below that present cosmologies are either consistent 
with the G/G found from the observations, or demand 
more QlQ and a negative w/cOrj^^. This situation 
appears to be fairly strictly either-or. In view 
of this, it is illustrative to view the history of 
geophysical papers requiring some nontidal earth 
acceleration, 

Munk & MacDonald (1960) and Dicke (1966), 

(1969) did not appeal to large nontidal acceler- 
ations in their theoretical analyses. The apparent 
cd\?,ervat?UmfXl values cried out for explanation, and 
several suggestions (see above; have been found 
by theoreticians. It may be unfair to suggest this, 
but the feeling persists that if it were suddenly 
possible to show aheervationallu that « 0> 

the theoretical analyses would soon follow to show 
that yes, after all, the earth-moon system tides do 
show apparent conservation of angular momentum. 
Against this we have several cosmologies, each of 
which demands at least as much apparent nontidal 
^/(jy (as G/G) as the data now provide. Tliese cosmolo- 
gies were not constructed in response to the non- 
tidal acceleration anomaly (as were the geophysical 
analyses) , and yet the need from all of them is for 
enough G/G to easily replace all the ih/avpj^q;. which 
the data can muster. This is too much of a common 
demand to dismiss as coincidence, and suggests that 
ft/G 0 is more likely than a significant di/fjjjjyj'j*. 


Primitive cosmologies, for all values of e, 
provide from figure 17; 

6/G = - 2.3 ± 1.5 X 10"”yr“‘ (65) 

Theories consistent with the figure 16 equations of 
condition Include: Hoyle &Narlikar (1974), Dicke 
(1957) & (1962) , Peebles & Dicke (1962), providing; 

H = 3st”^; G«T“^; G/G - - H (66) 

where T is the parametric age of the Universe, and 
H is the Hubble Constant. Sandage & Tammann (1975) 
obtain a Hubble Constant of: 

H - 55 ± 7 km/s/Mpe (67) 

which implies a fc/G from extrag^lactic observation: 

G/G - - 5.6 ± 0.7 X lO^^Vyr”^ (68) 

Equation (65) is over two standard deviations from 
the requirement of (68) and this is not good news 
for this class of cosmologies. On the other hand, 
the Brans & Dicke (1961) cosmology involves several 
parameters, and its implications probably bound &/G: 

0 > 6/G > -3H/(w4*2) > -5 X 10"^^yr“^ (69) 

With preferred values near -2 X 10”^^ yr"^ ; see 
Weinberg (1972) page 629. The result (65) is in 
adequate agreement with this theory. I am not 
competent to judge the roots of theoretical cosmol- 
ogy, but have been frequently advised that the 
entire class of primitive cosmologies is in deep 
trouble on other observational and theoretical 
considerations. No doubt there are those who 
disagree, so suffice it to say, the numerical 
results of this paper cast some observational doubt 
on the primitive cosmologies except Brans-Dicke. 

Dirac’s cosmology also demands that G/G - -H 
and (68), which is in good agreement with the 
solution from figure 17 line 2 on this cosmology; 

G/G = -5.1 ± 3.0 X 10~“yr“^ (70) 

The equations of condition implied by this cosmology 
are probably the best we have available, and this 
result is In good agreement with the finding of 
VanFlandem (1976), equations 4 and 5a of figure 16 
alone ; 

&/G = -5.8 ± 3.1 X 10"“ yr"’ (71) 

There are, in effect, two independent experimental 
ways to find 6/G represented in figure 16 equations 
of condition. One is given in (71) . Tlie other is 
to uj/a)Y^.jT= 0 and use equations 3a and 4 from 

figure 16. This uses the same smu'ce of as the 
reference, but an entirely independent source for 
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fl)/u) which i$ also independent of the moon. Solving 
these equations yields: 

f!/0 = - 6.9 ± 3.0 X 10"** yr"^ (72) 

It is difficult to ignore the patterns of 
strong agreement between various determinations of 

[1] from two independent lines of astronomical 
data from the solar system; [2] with the Hubble 
Constant determined from extragalactic data via the 
equations of the best available cosmology. [31 Tlie 
solution for G/G with the Hubble constraint (68) on 
ayiy current cosmology eliminates all significant 
nontidal earth acceleration [4] No geophys- 
ically assumed can eliminate the result (71) 

because VanFlandem's data is independent of the 
earth's rotation. 

There is only one interpretation which will 
satisfy all of the observations, and that is to 
simply accept the cosmological 6/G soloed from the 
data. If the alternative of a large nontidal 
earth acceleration is to be forced, hath of 

the following must hold. [a] tlie VanFlandern (1976) 
observation must be biased enough to obviate his 
result. His analysis is a very difficult under- 
taking, and he has had problems in eliminating 
systematic errors: compare VanFlandern (1975) and 
(1976). His latest results are very encouraging, 
and the solution uncertainties are now declining 
very rapidly because a full 21 year data span is 
breaking correlations resulting from the 18.6 year 
lunar nodal rate. His latest results appear very 
strong, and the numerical values have not changed 
significantly in the last year of analysis. [b] The 
tvue cosmology must require a zero or very small G/G 
compared with the Hubble rate. Every modern cosmol- 
ogy has established this relationship at an equal 
or greater level than that required to completely 
eliminate the observed as noted above. The 

Dirac cosmology may not be the final answer to every 
demand, but neither this nor any other specific 
cosmology needs to be adopted in order for this 
argument to follow. A cosmology intrinsically 
different from all of the recent suggestions must 
be advanced and justified before this constraint can 
be avoided. 

The acceptance of the cosmological &/G found 
here and in VanFlandern (1976) is consistent With 
all of the observations. To accept the alternative 
geophysical a3/tu.j.^rp it is necessary to abandon a 
basic observation set (Atomic Time lunar occulta- 
tions) , as well as the basic thrust of every major 
cosmology (the relation of G/G to the Hutble 
Constant) . This appears to be an unlikely, though 
perhaps not quite impossible, choice. Accepting 
the Hubble rate determination (68) as a further 
equation of condition added to figure 16 equations 
la - 5a probably provides the best available 
observational estimate of 

“/“TNT “ ^ ^ ^ ^ ^ 10"** yr"* 


The uncertainty of figure 17 line 2 is reduced by 
adopting the additional equation (68) instead of 
solving for all three parameters on the Dirac 
equations of condition and solar system observables 
only. Tlie larger uncertainty can be taken if the 
inclusion of extragalactic data is not desired, as 
can the solution of line 1 if the primitive cosmol- 
ogies are preferred. In any case, the i^^ 

consistent with zero, and there remains no need to 
find a geophysical explanation of any significant 
nontidal earth acceleration. 

It is concluded that there are strong reasons 
for accepting the cosmological ^J/G found from the 
observations, and consistent with the Hubble 
Constant. The alternative of a large nontidal 
acceleration of the earth's rotation is unlikely, 
and the earth-moon system tides appear to conserve 
angular momentum wholly within the tidal couple to 
about 3 X lO'^^yr"' or better. The tidal acceler- 
ations are observed to remain constant within 
uncertainties during the historical period. 

Summary of Principal Results 

• Tidal accelerations constant within standard error 

• Strong independent support for cosmological ft/G 

• Nontidal earth acceleration consistent with 0 

• Consistent with Dirac and Brans-Dicke cosmologies 
■ Ancient, modern, and extragalactic data agree 

Figure 18 

Principal Numerical Results 

ANCIENT OBSERVATION RESULTS 

- -30.0 ± 3.0 "cy’"^ Eclipses Only 

(ly/w - -24. 3 ± 2.0 X 10^^Tyr“^ Eclipses Only 

Ci/bi - -2 3.8 ± 2. 3 X 10"‘^^yr"'^ Timed Equinoxes 

D” - n - 1.61d)/w = +9.10 ± 0.3 


RESULTS FROM ALL ASTRONOMICAL DATA 

-27.2 ±1.7 *' cy*"^ 
b)fbi -ll\e ± 1.1 X 10^^^yr“^ 

AT = 20 + 114 - T + 35.55 - ±0 AT^ . 3T^ + 50 

ft == +4.39 "ey^^ ft - -.10 "cy“2 

6 / G = -2 . 3 ± 1 . 5 X 10”^^ yr"^ Primitive Cosmology 
G/G = -5.1 ± 3. O X 10““yr"V Dirac Cosmology 
Ui/tDxNT = +2 . 5 ± 2 . 5 TO +9 . 2 X 2 . 5 X 10' y r' ^ 
fSolve G/G± + G/G = 0 4 
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